Hi Everyone. As part of my project, I did pre-processing steps on my Radarsat-2 image which includes (Matrix Generation, Polarimetric Decomposition,Terrain correction). Then, I am extracting urban area using thresholding algorithm using python. Then, I have to show that area on Google Earth.I am posting my code here. Please let me know where I am going wrong.
import gdal
import matplotlib.pyplot as plt
import numpy as np
import osr
import cv2
gtif = gdal.Open(‘geosample1.tif’)
srcband = gtif.GetRasterBand(1)
trans = gtif.GetGeoTransform()
proj = gtif.GetProjection()
if proj != 0:
print(“1”)
else:
print(“Hi”)
arr = srcband.ReadAsArray()
print(“urban area”)
for i in range(arr.shape[0]):
for j in range(arr.shape[1]):
if 55.4525 <= arr[i,j] <= 94.2696:
arr[i,j] = 255
else:
arr[i,j] = 0
plt.imshow(arr,cmap = ‘gray’)
import gdal
import osr
cols = arr.shape[1]
rows = arr.shape[0]
driver = gdal.GetDriverByName(“GTiff”)
sr = osr.SpatialReference()
sr.ImportFromEPSG(4326)
sr_wkt = sr.ExportToWkt()
outRaster = driver.Create(‘E:\tanvi_iirs_work\Task 3\Images\output5file.tif’,cols,rows,1,gdal.GDT_Float32)
outRaster.SetGeoTransform(trans)
if outRaster:
print(‘Updating Projection’)
res = outRaster.SetProjection(sr_wkt)
if res != 0:
print(‘Setting Projection Failed’)
outband = outRaster.GetRasterBand(1)
outband.WriteArray(arr)
outband.FlushCache()
outRaster = None
else:
print(‘couldnt open’)