[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