[GRASS-user] mask in PERMANENT affect on other mapsets during reprojection

Jarekj jarekj at amu.edu.pl
Mon Jun 30 00:43:32 EDT 2008


So as I understand the problem will heppen when I have mask in mapset in 
source which is named as my target mapset in target location: ussually 
it will concern PERMANENT?

Thanks for information
Jarek

Glynn Clements pisze:
> [CC to grass-dev]
>
> Jarekj wrote:
>
>   
>> The situation is as follows:
>> In location GRS mapset PERMANET I have MASK (with projection of 2180 
>> based on GRS elipsoid)
>> in second mapset let say jarekjaI have map which is to be reprojected to 
>> other location (WGS). In that mapset no MASK is present.
>> In location WGS mapset PERMANENT  I reporject map from location GRS and 
>> I have only fragment allowed by MASK from GRS PERMANENT It is intended 
>> or some bug?
>>     
>
> The MASK file has to exist in the current mapset in the current
> location.
>
> However: r.proj keeps switching the current location between the two
> locations. It switches to the source location when opening and reading
> the input map, and switches back to the target location when opening
> and writing the output map.
>
> Where it looks for the mask is determined by which is the current one
> at the point that it looks. This point isn't well defined; it normally
> happens as a side effect of some other function.
>
> AFAICT, it uses the source location when looking for the mask. 
> However, it doesn't change the mapset; i.e. it uses the current mapset
> (not the mapset= value), which may not even exist in the source
> location.
>
> This is quite clearly a bug; it should switch mapsets when it switches
> locations. At present the code only uses the mapset= value when
> opening the input map, overlooking the fact that it is also used
> internally, e.g. when looking for the mask.
>
>   


More information about the grass-user mailing list