use "gdal_translate" function in GDAL
# get the corner xy coordinates
gdalinfo "filename"
# define corner xy coordinates
gdal_translate -gcp 0.0 0.0 488730.762621 3836122.524240 -gcp 7132.0 0.0 558280.869646 3820108.908850 -gcp 0.0 7000.0 473142.808142 3767848.968870 -gcp 7132.0 7000.0 542618.864253 3751940.046870 -of "GTiff" [input filename] [temporary filename]
--memo
(-gcp [former x coordinate] [former y coordinate] [defined x coordinate] [defined y coordinate])
(upper left coordinate, upper right coordinate, lower left coordinate, lower right coordinate)
--memo
# define EPSG and make it into Geotiff
gdalwarp -t_srs "EPSG:32652" [temporary filename] [output filename]
