[GRASS-user] Copy raster map from one unprojected location to another

Nikos Alexandris nik at nikosalexandris.net
Thu Mar 19 05:11:51 PDT 2015


Nikos Alexandris wrote:

..

> So, why go into checking anything regarding PROJ_INFO or PROJ_UNITS
> if the "-o" flag is instructed?
> Completely untested, but maybe the script needs a modification (?), 
> as follows:
>
>     # 
> -----------------------------------------------------------------------
>     if flags['o']:
>         grass.warning(_("Overriding projection check (using current
> location's projection)!"))
>
>     else:
>
>         diff_result_1 = diff_result_2 = None
>
>         proj_info_file_1 = 'PROJ_INFO'
>         proj_info_file_2 = os.path.join(mset_dir, '..', 'PERMANENT',
> 'PROJ_INFO')
>
>         if not 
> grass.compare_key_value_text_files(filename_a=proj_info_file_1,
>                                                   
> filename_b=proj_info_file_2,
>                                                   proj=True):
>             diff_result_1 = diff_files(proj_info_file_1, 
> proj_info_file_2)
>
>         proj_units_file_1 = 'PROJ_UNITS'
>         proj_units_file_2 = os.path.join(mset_dir, '..', 'PERMANENT',
> 'PROJ_UNITS')
>
>         if not 
> grass.compare_key_value_text_files(filename_a=proj_units_file_1,
>                                                   
> filename_b=proj_units_file_2,
>                                                   units=True):
>             diff_result_2 = diff_files(proj_units_file_1, 
> proj_units_file_2)
>
>         if diff_result_1 or diff_result_2:
>
>             if diff_result_1:
>                 grass.warning(_("Difference between PROJ_INFO file of
> packed map "
>                                 "and of current
> location:\n{diff}").format(diff=''.join(diff_result_1)))
>             if diff_result_2:
>                 grass.warning(_("Difference between PROJ_UNITS file
> of packed map "
>                                 "and of current
> location:\n{diff}").format(diff=''.join(diff_result_2)))
>             grass.fatal(_("Projection information does not match. 
> Aborting."))
>     # 
> -----------------------------------------------------------------------

A quick test after the modification works for me.  Though, not properly 
tested... :-)

Nikos


More information about the grass-user mailing list