plotting OGR_GMT polygons with psxy
|Target version:||Candidate for next bugfix release|
|Affected version:||5.x-svn||Platform:||Mac OS X|
Dear GMT gurus,
I have noticed an error in newer versions of GMT > 5.1.x when attempting to plot OGR_GMT format polygon files. As a test of the different GMT versions, I made a very simple plotting script:
#!/bin/bash output=africa_vegmap.ps # ----------------------------- psbasemap -R-25/-34/55/25r -JA20/0/20 -D-25/-34/55/25r -B0 -P -K > $output psxy potveg_biomes.gmt -R -J -A -fg -aG=gmtrgb -L -O -P >> $output ps2raster -A -Tf $output rm $output
Both the correct (with GMT 5.1.2) and incorrect (with GMT 5.3.0_r17024) versions of the output are attached to this bug report, along with the OGR_GMT polygon file I am using for plotting. Basically, the older version of GMT recognizes the OGR_GMT polygon file correctly and plots each polygon in the correct color. The newer version of GMT somehow thinks that all of the polygons should be painted the same color.
I recompiled both versions of GMT today, so they are linked to exactly the same libraries, including, e.g., GDAL.
Thank you very much for your help resolving this!
- Status changed from In Progress to Resolved
- % Done changed from 0 to 100
In a fix to deal with a related problem we ended failing to parse the @P and @H properly so that once a hole was detected everything else was flagged as holes; hence everything became blue after the last hole. Fixed in r17150 - give it a try.