Messages

The TIGER data contains not only shapefile sets (edges.dbf, edges.prj, edges.shp, edges.shp.xml, edges.shx), but also database only file sets (featnames.dbf, featnames.dbf.xml). If I understand this correctly the featnames.dbf contains identification information for records contained in edges.shp and they are associated by the TLID. I know how to load the edges.shp file into PostgreSQL/PostGIS (shp2pgsql, psql) and I think I can load the featnames.dbf (shp2pgsql, psql). Now how do I associate the TLID so that I can generate shapefiles that contain selected additional data from the featnames table.

Thanks

Just FYI. I recently used ver 1.0.21 for the transportation files and discovered that files are now named 2009 while topo process still expects 2008. A simple work around is to download the files, and manually rename before continuing.

Also if you only have a single contour file, don't check "Combine..." as everything gets deleted except for 2 .prj files.

Quite a handy tool though.

Yeah, I was fighting with the very same issue. After upgrading everything and no success, I turned to the web. Looks like you posted around the time I started upgrading things.  Moved my maps to a shorter directory path (C:\maps) and all is better now.

Is it possible to find out what the erroneous projection is and then correct for it?   I thought it was just an incorrectly specified datum, (Guam  1963) but after converting it still seems that it is still in error.  The error appears to be different as latitude increases.

I noticed that you have recently uploaded the map of the Commonwealth of the Northern Marianas Islands TOPO Map.  Thank you very much. It looks good.

In a post somewhere on this forum, that I can no longer find, I recall you saying that you were working on a map of Guam. The thread was in reference to US Territories and Protectorates raster files. Are you planning to produce such a map? I am very much interested.

I attempted this myself about a year ago and ended up with a less than satisfactory product. Since some time has passed, I assumed that data sources would be better and newer tools will have simplfied the process. I'm disappointed to say that I'm having a more difficult time now than I originally did, but I'm not giving up just yet.

Roy

Dan,
Thank you so much for the program update.  Processing the Guam data results in an error "Run-time error '70': Permission denied".  Arizona works fine. I wonder if it is due the special characters in some of the names such as "Hagåtña".

One other problem I still have is that I need to change positions of the gpx points. Or can I project the shp file?

Roy

I'm not making much progress at all.  I was hoping to figure some of this out myself, but now I'm back asking for help again.  I got the GNIS file, converted it to gpx, using the "USGS GNIS to GPX Converter", examined the resulting gpx file, and discovered the elevations were in meters.  Converted the meters to feet, and then added to my .mp  contours only to discover that the points were all shifted again (datum problems).  Translated the gpx from Guam_1963 to NAD83.  They are now cottectly positioned.  At this point I am stuck again, I cannot figure out how to use gpsmapedit to generate shp files from the gpx as instructed above.  I thought the gpx should be loaded into the postgis database and then be extracted as a shp file, but I am unable to figure out how to do this.  Any help would be greatly appreciated.
Roy

Dan,
I've annotated the lines on the attached jpeg.  I guess my color choices didn't provide enough contrast.
Roy

Thank you for the GNIS file.  I looks like I now could figure out how to do other areas from the original zips with a little perl (how did you greate this file?).

Oz,
I assume this is the first step or two that GNIS Process does.  Per part 5 of the tutorial GNIS Process will make a shp file.  What steps do I need to do to turn this csv into a shp?  I'm guessing it needs to get put into the postgis database and then pulled out as a shp file using pgsql2shp, but I'm at a loss with how to do it.

Roy

Dan,
Thank you for your efforts.

This is still somewhat puzzling to me but I have an idea. When importing the converted edges file you provided (census) seems to be about 50 feet different from the data I got from the NHD.

As for the contour data, I finally have figured out how to use FWTools (I think, at least it produces output).  Using gdalwarp, I assumed the geotiff data was in Guam_1963 datum and converted it to NAD83 (even though the seamless server claims to be NAD83). Then generated contours with dem2topo.  Loaded the mp file and both the census edges and NHD data into GPSMapEdit (census is in black and NHD is in blue).  The HND data seems to match the contours quite well.

I have two data points from a Garmin 60csx that seem to match quite favorably with the contour data.  I am still trying to get some better points to verify the contour data.

What I still don't understand is why the edges from census and flowlines from NHD are so different.

Can you sed some light on my use of gdalwarp?  I used:
gdalwarp -s_srs EPSG:4675 -t_srs NAD83  input.tif  output.tif

Thanks,
Roy

Oz,
It appears that you are the author of GNIS Process.

I am trying to build a map of Guam (territory of the US) and maybe some of the other pacific islands.  The current GNIS Process only lists the 50 states, but I know data for the US territories exists at geonames.usgs.gov.  Could these be added to GNIS Process or could you share the process that GNIS Process does?  I'm thinking I could do it by hand.

BTW When running Hawaii as a test,  I get a "Run-time error '62': Input past end of the" in a dialog while the status says "Extracting GNIS Zip File".  Testing with Minnesota runs fine.

Thanks,
Roy

Here's the Dr Watson log from the crash.  Made a new one with less stuff running in Windows.

I've been fighting with this on multiple fronts.  Hopefully I'm keeping it all straight.
The NHD data link is:
ftp://nhdftp.usgs.gov/NHD136010.zip

I have loaded the NHDFlowline.shp fie into GPSMapEdit specifying NAD83 as Datum

then from TIGER 2007 data:
ftp://ftp2.census.gov/geo/tiger/TIGER2007FE/66_GUAM/66010_Guam/
I imported fe_2007_66010_edges.shp specifying Guam 1963 as Datum

Both of these shp files seem to compare well (within about 100 ft), but are significantly different than the contour data.

Will look into finding a reference point to try and determine which is correct.  Is it possible to tell dem2topo that the tif is in datum other than WGS84?

I realize that Fenwick started this thread, but since I am having what seems to be the same problem, is it ok for me to post here also?  I don't want to violate any forum etiquette.  If this is not acceptable let me know and I'll start another thread.

My error info:
processing guam NHD from TIGER 2007 data
ftp://ftp2.census.gov/geo/tiger/TIGER2007FE/66_GUAM/66010_Guam/

Could not find the dump from pgsql2shp, but did locate a log in the ...\PostgreSQL\8.3\data\pg_log directory
2008-10-18 11:45:40 PDT LOG:  database system was shut down at 2008-10-18 11:45:22 PDT
2008-10-18 11:45:40 PDT LOG:  loaded library "$libdir/plugins/plugin_debugger.dll"
2008-10-18 11:45:41 PDT FATAL:  the database system is starting up
2008-10-18 11:45:41 PDT LOG:  database system is ready to accept connections
2008-10-18 11:45:41 PDT LOG:  autovacuum launcher started
2008-10-18 11:45:42 PDT LOG:  loaded library "$libdir/plugins/plugin_debugger.dll"
2008-10-18 11:47:25 PDT LOG:  loaded library "$libdir/plugins/plugin_debugger.dll"
2008-10-18 11:47:25 PDT NOTICE:  CREATE TABLE will create implicit sequence "nhd_area_gid_seq" for serial column "nhd_area.gid"
2008-10-18 11:47:26 PDT NOTICE:  CREATE TABLE / PRIMARY KEY will create implicit index "nhd_area_pkey" for table "nhd_area"
2008-10-18 11:47:35 PDT LOG:  loaded library "$libdir/plugins/plugin_debugger.dll"
2008-10-18 11:47:45 PDT LOG:  loaded library "$libdir/plugins/plugin_debugger.dll"
2008-10-18 11:48:02 PDT LOG:  could not receive data from client: No connection could be made because the target machine actively refused it.

   
2008-10-18 11:48:02 PDT LOG:  unexpected EOF on client connection

I don't think this is what we need. Will try to figure out why there is no error file even though windows says it's creating one.

Funny,
I get exactly the same thing (you beat me to the post).  Only difference is I am running Win2K natively.

I installed PostgreSQL 8.3 and PostGIS 1.3.3 for PostgreSQL 8.3

I'm guessing there are fields or values from the original shp files that pgsql2shp is not set up to handle.

Maybe if we can get some documentation on pgsql2shp we could help trouble shoot this.

Updated...
I have found some very good documentation on PostGIS.  Guess I should have looked around first.  No guarantees I can figure this out, but I'm going to give it a try.