[GRASSLIST:5444] problems with r.in.srtm

samuel cavalcante rioapaporis at yahoo.com.br
Wed Jan 19 14:44:47 EST 2005


hi,
i tried to import .hgt files using r.in.srtm but it didn´t work...r.in.srtm instead produces 3 files ... then i did run r.in.gdal to import the  .bil file, but it didn´t recognize well the projection...i have another rasters i want to overlay but even if they (apparently) are in the same projection (lat/long wgs84) they don´t match...
can i change the projection of a raster without changing the location?
thanks
samuel
 

Carlos Henrique Grohmann <guano at usp.br> wrote:
I'm still with problems with my vect maps, so now I'll post all msgs.

I have this USGS Open File Report (97-470D, 1999) with a shape file of the
political divisions of South America. I want to import it and make a vect map
of South America without divisions. So I've imported it like this:


>v.in.ogr dsn=/mnt/cdrw_dvdrom/SAM/VIEWS/GEO output=cont5 layer=POL6_G
min_area=0.0001 >type=boundary snap=-1 -o 
>Over-riding projection check.
>Proceeding with import...
>Layer: POL6_G
>-----------------------------------------------------
>Building topology ...
>Registering lines: 
>337 primitives registered 
>Building areas: 
>
>196 areas built 
>177 isles built
>Attaching islands: 
>Attaching centroids: 
>Topology was built.
>Number of nodes : 319
>Number of primitives: 337
>Number of points : 0
>Number of lines : 0
>Number of boundaries: 337
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 177
>Number of incorrect boundaries : 93
>Number of areas without centroid : 196
>-----------------------------------------------------
>
*** stuff ***
>
>Building topology ...
>Building areas: 
>
>196 areas built 
>182 isles built
>Attaching islands: 
>Topology was built.
>Number of nodes : 319
>Number of primitives: 341
>Number of points : 0
>Number of lines : 4
>Number of boundaries: 337
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 182
>Number of areas without centroid : 196
>
>Layer: POL6_G
>-----------------------------------------------------
>-----------------------------------------------------
>Building topology ...
>Registering lines: 
>337 primitives registered 
>Building areas: 
>
>196 areas built 
>182 isles built
>Attaching islands: 
>

And I got my map, with all countries in South America. Then I try to dissolve
all commom boundaries with:


>v.extract input=cont5 output=cont6 type=boundary layer=1 new=1 -d 
>Building topology ...
>Registering lines: 
>333 primitives registered 
>Building areas: 
>
>196 areas built 
>182 isles built
>Attaching islands: 
>Attaching centroids: 
>Topology was built.
>Number of nodes : 319
>Number of primitives: 333
>Number of points : 0
>Number of lines : 0
>Number of boundaries: 333
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 182
>Number of areas without centroid : 196
>
>Writing attributes ...
>Removing duplicate centroids ...Building topology ...
>Registering lines: 
>
>333 primitives registered 
>Building areas: 


And I got the same map! No dissolved commom boundaries..

So I did:



>v.extract input=cont5 output=cont6 type=boundary layer=1 new=1 -d 
>Building topology ...
>Registering lines: 
>333 primitives registered 
>Building areas: 
>
>196 areas built 
>182 isles built
>Attaching islands: 
>Attaching centroids: 
>Topology was built.
>Number of nodes : 319
>Number of primitives: 333
>Number of points : 0
>Number of lines : 0
>Number of boundaries: 333
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 182
>Number of areas without centroid : 196
>
>Writing attributes ...
>Removing duplicate centroids ...Building topology ...
>Registering lines: 
>
>333 primitives registered 
>Building areas: 

All the same, the boundaries are still there.


Once I get an error like this:

>bgerror failed to handle background error.
> Original error: too many nested calls to Tcl_EvalObj (infinite loop?)
> Error in bgerror: too many nested calls to Tcl_Eval (infinite loop?)

and then GRASS crashed.

Pleas help. I don't know what to do.

I'll leave the shp file in 

http://www.igc.usp.br/pessoais/guano/sam.html


thanks


-- 
+----------------------------------------------------------+
Carlos Henrique Grohmann - Guano 
Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - guano at usp dot br
+----------------------------------------------------------+


		
---------------------------------
Yahoo! Acesso Grátis - Internet rápida e grátis. Instale o discador do Yahoo! agora.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20050119/996542e3/attachment.html


More information about the grass-user mailing list