R: Re: R: Re: [GRASS-user] Input raster map is outside current region

[Micha Silver]

On 03/11/2010 04:49 PM, ciccpaol at libero.it wrote:
> Hi Micha,
>
> i fixed my previous problem...i can see the map on the display (it was because
> of the wrong resolution setting: in fact i had to add to the r.proj command the
> parameter res=<resolution of the raster converted into the map units of the
> target location>)
>    
ׂGood. You probably want to set the res even before r.proj, with g.region 
res=xxx
> I can see my raster map, but the raster changed shape: all the printed squares
> of the utm system became rectangles....mmm
>    
That's what reprojection does!
> i guess it has got smthing to do with the resolution setting...
>
>
>    
>> ----Messaggio originale----
>> Da: micha at arava.co.il
>> Data: 11/03/2010 14.47
>> A: "ciccpaol at libero.it"<ciccpaol at libero.it>
>> Cc:<grass-user at lists.osgeo.org>
>> Ogg: Re: R: Re: [GRASS-user] Input raster map is outside current region
>>
>> On 03/11/2010 03:34 PM, ciccpaol at libero.it wrote:
>>      
>>> Hi ,
>>>
>>> thanks Micha for your help! I tried the steps you suggested me, but they
>>>        
> didnt
>    
>>> work properly; after trying that, i found an article about that
>>> http://wiki.gfoss.it/index.php/Riproiezioni_di_dati_raster
>>> that is more or less the same you suggested me.
>>>
>>> Working those steps, I finally got the right region extent. In fact when
>>>        
> I'm
>    
>>> inside the target location and I'm displaying the reprojected raster map,
>>>        
> the
>    
>>> coordinates values shown on the bottom of the page are right, but I can't
>>>        
> see
>    
>>> the map! The page is all blank.
>>>
>>>
>>>        
>> I guess it's time to ask for more information: Can you post the outputs
>> for *both locations* of:
>>
>> (Source location)
>> g.region -p
>> v.info -g<src vector region as below>
>> r.info -g<original raster>
>>
>> (Target location)
>> g.region -p
>> v.info -g<projected vector region>
>> r.info -g<final projected raster which doesn't appear>
>>
>>      
>>>> What you do here is:
>>>>
>>>> - go back to the source location and run
>>>> v.in.region out=src_region
>>>> This creates a vector polygon of the bounds of the current region.
>>>>
>>>> - now switch to the new, target region and do
>>>> v.proj location=<src_location>   mapset=<src_mapset>   in=src_region
>>>> out=tar_region
>>>> to reproject that rectangular polygon to the target region . Next run
>>>> g.region vect=tar_region
>>>> And the region in your target location should match where the raster
>>>> wants to be.
>>>>
>>>> - Now try your r.proj comand
>>>>
>>>>          
>>>
>>>        
>>>>>
>>>>>            
>>>> -- 
>>>> Micha Silver
>>>> http://www.surfaces.co.il/
>>>> Arava Development Co.  +972-52-3665918
>>>>
>>>>
>>>>
>>>>          
>>> _______________________________________________
>>> grass-user mailing list
>>> grass-user at lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>>
>>> This mail was received via Mail-SeCure System.
>>>
>>>
>>>
>>>        
>>
>> -- 
>> Micha Silver
>> http://www.surfaces.co.il/
>> Arava Development Co.  +972-52-3665918
>>
>>
>>      
>
>
> This mail was received via Mail-SeCure System.
>
>
>    


-- 
Micha Silver
http://www.surfaces.co.il/
Arava Development Co.  +972-52-3665918