|
From: Arlindo da S. <da...@al...> - 2007-11-26 02:03:39
|
Graziano,
I have an initial implementation of your shapefile "draw polyf"
implemented as a GrADS extension. Just check out the Grads module
% gacvs co Grads
% cd Grads/extensions/shape
% make
% perldoc shape
% setenv GAUDXT ./shape.udxt
For convenience I have incorporated libshp with the sources, at least
for now. Look at file utShape.gs for an example. You will need the
gshhs_land.??? files for running this test, and a GrADS build
supporting UDXTs. You should be able to build the sources you just
checked out.
A few comments:
1) I rearranged your code a bit, grouping together portions that could be reused
2) I added "shp_lines" for my own benefit, it may not be needed in the end as
polyf should handle non-closed polygons as well. At this point,
though, they do not produce the same output.
3) I am not sure whether "shp_polyf" works as you intended. The
polygon fills look funny.
4) The current sources let GrADS handle clipping. Your old function
NormalLon (renamed here StdLon) and the InsideDomain() functions are
effectively disabled. You may have good reasons for re-enabling them.
5) The documentation needs to be filled in.
6) It would be nice to have a nice set of simple examples, with small files.
7) It appears to me that one should be able to write shapefiles with
the contents of
the GrADS "metafile" (with coordinates expressed in lat/lon which I
believed can be accomplished with the appropriate conversion). I went
ahead and put in the skeleton from my "gxyat" extension which grabs
the metafile buffer from memory and translates it. This work would be
better completed by someone with shapefile familiarity. Better yet,
with means of testing it. I have no experience (or personal interest
for that matter) in GIS software, but I realize it would be important
to reach out to the GIS community. Any volunteers to take this on? (I
can help with the metafile translation if needed.)
Efren: have you been able to make progress in this area?
Please let me know if you have any comments.
Arlindo
--
Arlindo da Silva
da...@al...
|