[Gdal-dev] tigerpoly.py - how to use it?

[Jack Twilley]

I am an extreme newbie to GDAL, so please be gentle.

I discovered pymod/samples/tigerpoly.py and tried to use it.

I downloaded the TIGER/Line 2002 data for Contra Costa County in
California and have the following files:

total 34580
drwxr-xr-x    2 jmt  jmt       512 Nov 23 23:39 .
drwxr-xr-x  103 jmt  jmt      8192 Nov 24 10:55 ..
-rw-r--r--    1 jmt  jmt  12345480 Nov 20  2002 TGR06013.RT1
-rw-r--r--    1 jmt  jmt   3426150 Nov 20  2002 TGR06013.RT2
-rw-r--r--    1 jmt  jmt    101700 Nov 20  2002 TGR06013.RT4
-rw-r--r--    1 jmt  jmt    777954 Nov 20  2002 TGR06013.RT5
-rw-r--r--    1 jmt  jmt    228930 Nov 20  2002 TGR06013.RT6
-rw-r--r--    1 jmt  jmt     37240 Nov 20  2002 TGR06013.RT7
-rw-r--r--    1 jmt  jmt     14668 Nov 20  2002 TGR06013.RT8
-rw-r--r--    1 jmt  jmt   3030116 Nov 20  2002 TGR06013.RTA
-rw-r--r--    1 jmt  jmt     30256 Nov 20  2002 TGR06013.RTC
-rw-r--r--    1 jmt  jmt   2050240 Nov 20  2002 TGR06013.RTH
-rw-r--r--    1 jmt  jmt   6924204 Nov 20  2002 TGR06013.RTI
-rw-r--r--    1 jmt  jmt    671771 Nov 20  2002 TGR06013.RTP
-rw-r--r--    1 jmt  jmt       312 Nov 20  2002 TGR06013.RTR
-rw-r--r--    1 jmt  jmt   2429810 Nov 20  2002 TGR06013.RTS
-rw-r--r--    1 jmt  jmt   1932903 Nov 20  2002 TGR06013.RTT
-rw-r--r--    1 jmt  jmt   1117088 Nov 20  2002 TGR06013.RTZ

When I run the tigerpoly.py script using input and output directories
for arguments, I get the following shapefile:

-rw-r--r--    1 jmt  jmt  4962200 Nov 23 23:41 out.dbf
-rw-r--r--    1 jmt  jmt  4475352 Nov 23 23:41 out.shp
-rw-r--r--    1 jmt  jmt   114444 Nov 23 23:41 out.shx

shpdump tells me that this shapefile is of type polygon and contains
14293 shapes.  When I examined the output with the rendering tool I
normally use for shapefiles, it appears to have all the possible
layers jammed together.  I see water polygons, land boundaries, major
roads, minor roads, water features, everything all together.

Is there a way to generate individual shapefiles for each layer
instead of one shapefile containing everything?  I read through the
script and didn't see any way to set options to only output a single
layer per file or anything.  If I'm missing something. please let me
know as this would be terribly useful.

Jack.
-- 
Jack Twilley
jmt at twilley dot org
http colon slash slash www dot twilley dot org slash tilde jmt slash
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 186 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/gdal-dev/attachments/20031124/407c9c94/attachment.bin