[gradsusr] [EXTERNAL] Re: Plotting shapefile elements

Mon Mar 7 13:34:57 EST 2022

There is an old demo scripts in the library:
ftp://cola.gmu.edu/grads/scripts/shp_demo1.gs (you might have to use wget to download this – some browsers are not supporting ftp urls anymore)

That shows how you would query the shapefile database and use that information to draw individual shapes and colorize them accordingly. Should be enough to get you started.
--Jennifer

--
Jennifer Miletta Adams
Senior Scientific Software Developer
Goddard Earth Sciences Data and Information Services Center (GES DISC)
NASA/GSFC, Code 610.2
Building 32, Room S159

From: gradsusr <gradsusr-bounces at gradsusr.org> on behalf of Jeff Duda <jeffduda319 at gmail.com>
Reply-To: GrADS Users Forum <gradsusr at gradsusr.org>
Date: Monday, March 7, 2022 at 1:02 PM
To: GrADS Users Forum <gradsusr at gradsusr.org>
Subject: [EXTERNAL] Re: [gradsusr] Plotting shapefile elements

Consult the GrADS user's page on shapefiles first: http://cola.gmu.edu/grads/gadoc/shapefiles.html<https://gcc02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fcola.gmu.edu%2Fgrads%2Fgadoc%2Fshapefiles.html&data=04%7C01%7Cjennifer.m.adams%40nasa.gov%7Cdf9b8eed8f5542b4d26008da0064930e%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C637822729215027072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=wVrwmO3ec4DhAhzM4qrpfAsMB9qPeV%2FkS1tLqzRujD4%3D&reserved=0>

I recommend writing each set of product types to a separate shapefile using 'gxout shp' and 'shpattr' to control which elements are draw to which file. Then plot all elements of each shapefile separately using traditional means.

Jeff Duda

On Mon, Mar 7, 2022 at 4:16 AM Jason M. Warren <jason at solargrafix.com<mailto:jason at solargrafix.com>> wrote:
Hello,

I am trying to plot a shapefile which contains multiple elements.  I
queried the shp and dbf files as follows:

ga-> q shp current_warnings
Shapefile Type=Polygon #Shapes=4 XBounds=-88.9:-86.86 YBounds=29.38:37.16
0:  Polygon  parts=1  vertices=5  XBounds=-88.05:-86.86
YBounds=29.38:30.01  ZBounds=0:0  MBounds=0:0
1:  Polygon  parts=1  vertices=19  XBounds=-88.49:-87.05
YBounds=36.5:37.16  ZBounds=0:0  MBounds=0:0
2:  Polygon  parts=1  vertices=8  XBounds=-88.9:-88.36
YBounds=31.43:31.63  ZBounds=0:0  MBounds=0:0
3:  Polygon  parts=1  vertices=11  XBounds=-88.35:-88.01
YBounds=32.05:32.25  ZBounds=0:0  MBounds=0:0

ga-> q dbf current_warnings
RECORD#,PHENOM,SIG,WARNID,WFO,EVENT,ISSUANCE,EXPIRATION,URL,MSG_TYPE,PROD_TYPE

0,MA,W,O.NEW.KMOB.MA.W.0008.0109022
         2070225.670784,KMOB,0008,2022-01-09 19:02:25,2022-01-09
21:00:00,http://forecast.weather.gov/product.php?site=MOB&issuedby=MOB&product=SMW,SMW,Marine<https://gcc02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fforecast.weather.gov%2Fproduct.php%3Fsite%3DMOB%26issuedby%3DMOB%26product%3DSMW%2CSMW%2CMarine&data=04%7C01%7Cjennifer.m.adams%40nasa.gov%7Cdf9b8eed8f5542b4d26008da0064930e%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C637822729215027072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=qtkuA5rPRHeOdstj5pgdMeh9b91qSIwD7julI4WBZlE%3D&reserved=0>
Warning
1,FF,W,O.CON.KPAH.FF.W.0005.0109022
         2055627.670773,KPAH,0005,2022-01-09 17:56:27,2022-01-09
21:15:00,http://forecast.weather.gov/product.php?site=PAH&issuedby=PAH&product=FFS,FFS,Flash<https://gcc02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fforecast.weather.gov%2Fproduct.php%3Fsite%3DPAH%26issuedby%3DPAH%26product%3DFFS%2CFFS%2CFlash&data=04%7C01%7Cjennifer.m.adams%40nasa.gov%7Cdf9b8eed8f5542b4d26008da0064930e%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C637822729215027072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=0wI1U7RhZ%2B%2FD5wHYxsQYkay1hmS2Qkf8NWrw7DS9emg%3D&reserved=0>
Flood Warning
2,SV,W,O.NEW.KMOB.SV.W.0003.0109022
         2072922.670791,KMOB,0003,2022-01-09 19:29:22,2022-01-09
20:15:00,http://forecast.weather.gov/product.php?site=MOB&issuedby=MOB&product=SVR,SVR,Severe<https://gcc02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fforecast.weather.gov%2Fproduct.php%3Fsite%3DMOB%26issuedby%3DMOB%26product%3DSVR%2CSVR%2CSevere&data=04%7C01%7Cjennifer.m.adams%40nasa.gov%7Cdf9b8eed8f5542b4d26008da0064930e%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C637822729215027072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=uMifRxa39KLdSa%2BYeLHtWWhEjKNGQByDDZCfLZJdq0E%3D&reserved=0>
Thunderstorm Warning
3,SV,W,O.CON.KMOB.SV.W.0001.0109022
         2072321.670789,KMOB,0001,2022-01-09 19:23:21,2022-01-09
19:45:00,http://forecast.weather.gov/product.php?site=MOB&issuedby=MOB&product=SVS,SVS,Severe<https://gcc02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fforecast.weather.gov%2Fproduct.php%3Fsite%3DMOB%26issuedby%3DMOB%26product%3DSVS%2CSVS%2CSevere&data=04%7C01%7Cjennifer.m.adams%40nasa.gov%7Cdf9b8eed8f5542b4d26008da0064930e%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C637822729215027072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=X%2FmetzjrgDAYdNAZSt9cJg6%2FPipzLg3uW1x0igy%2FE5c%3D&reserved=0>
Thunderstorm Warning

However, I can not figure out how to plot this with different colors for
each element (SVR, SMW, etc...).  I have tried the following code which
plots everything in the shapefile...

'set rgb 1  255 0 0 255'
'set shpopts -1'
'set line 1 1 1'
'draw shp current_warnings'

But I need to know how to draw each kind of element separately (SVR,
SMW, etc...).

Thanks,

Jason Warren

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr<https://gcc02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr&data=04%7C01%7Cjennifer.m.adams%40nasa.gov%7Cdf9b8eed8f5542b4d26008da0064930e%7C7005d45845be48ae8140d43da96dd17b%7C0%7C0%7C637822729215027072%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=i%2BUCssLExUr%2FM3JLiZKOkWoMlxIN3eT8P7bdBBqQZ20%3D&reserved=0>

--
Jeff Duda, Research Scientist
University of Colorado Boulder
Cooperative Institute for Research in Environmental Sciences
NOAA/OAR/ESRL/Global Systems Laboratory
Boulder, CO

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20220307/bd3ffb83/attachment-0001.html>