增加军民港口文件区分
parent
a9f390d2f8
commit
a66f9ac220
|
|
@ -16,8 +16,6 @@ spacetySliceEnvPathExEPath=r"d:\ProgramData\anaconda3\envs\spacetySliceEnv\pytho
|
||||||
logPath=r"R:\TYSAR-德清院\A-预处理-未标注\A0-算法版本\AA\SpacetySliceDataTools\log\process.log"
|
logPath=r"R:\TYSAR-德清院\A-预处理-未标注\A0-算法版本\AA\SpacetySliceDataTools\log\process.log"
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
def existOrCreate(dirpath):
|
def existOrCreate(dirpath):
|
||||||
if not os.path.exists(dirpath):
|
if not os.path.exists(dirpath):
|
||||||
os.makedirs(dirpath)
|
os.makedirs(dirpath)
|
||||||
|
|
|
||||||
|
|
@ -17,6 +17,11 @@ portsliceSize=5000
|
||||||
shipsliceSize=1024
|
shipsliceSize=1024
|
||||||
BlockOverLayer=0.25
|
BlockOverLayer=0.25
|
||||||
|
|
||||||
|
def existOrCreate(dirpath):
|
||||||
|
if not os.path.exists(dirpath):
|
||||||
|
os.makedirs(dirpath)
|
||||||
|
|
||||||
|
|
||||||
def get_filename_without_ext(path):
|
def get_filename_without_ext(path):
|
||||||
base_name = os.path.basename(path)
|
base_name = os.path.basename(path)
|
||||||
if '.' not in base_name or base_name.startswith('.'):
|
if '.' not in base_name or base_name.startswith('.'):
|
||||||
|
|
@ -212,20 +217,6 @@ def stretchProcess(infilepath,outfilepath,strechmethod):
|
||||||
print("图像拉伸处理结束")
|
print("图像拉伸处理结束")
|
||||||
|
|
||||||
|
|
||||||
#切片模式
|
|
||||||
def getSlicePoints(h):
|
|
||||||
n = int(math.floor((h - 1024) * 1.2 / sliceSize))
|
|
||||||
step=int(math.ceil((h-1024)/n))
|
|
||||||
ti=list(range(0,h-1024,step))
|
|
||||||
ti.append(h-1024)
|
|
||||||
# 评价重叠率
|
|
||||||
movelayer=[]
|
|
||||||
for i in range(len(ti)-1):
|
|
||||||
movelayer.append((ti[i] + 1024 - ti[i + 1]) / 1024 * 100.0)
|
|
||||||
print("重叠率:",movelayer)
|
|
||||||
return ti
|
|
||||||
|
|
||||||
|
|
||||||
def getsliceGeotrans(GeoTransform,Xpixel,Ypixel):
|
def getsliceGeotrans(GeoTransform,Xpixel,Ypixel):
|
||||||
XGeo = GeoTransform[0]+GeoTransform[1]*Xpixel+GeoTransform[2]*Ypixel
|
XGeo = GeoTransform[0]+GeoTransform[1]*Xpixel+GeoTransform[2]*Ypixel
|
||||||
YGeo = GeoTransform[3]+GeoTransform[4]*Xpixel+GeoTransform[5]*Ypixel
|
YGeo = GeoTransform[3]+GeoTransform[4]*Xpixel+GeoTransform[5]*Ypixel
|
||||||
|
|
@ -254,9 +245,9 @@ def getNextSliceNumber(n,sliceSize,overlap=0.25):
|
||||||
return newN,ti
|
return newN,ti
|
||||||
|
|
||||||
def sliceShipDataset(rootname,im_data,src_im_data, im_Geotrans, im_proj, outfolder):
|
def sliceShipDataset(rootname,im_data,src_im_data, im_Geotrans, im_proj, outfolder):
|
||||||
binfolder=os.path.join(outfolder,"unit8binfolder")
|
binfolder=os.path.join(outfolder,"舰船","unit8binfolder")
|
||||||
pngfolder=os.path.join(outfolder,"pngfolder")
|
pngfolder=os.path.join(outfolder,"舰船","pngfolder")
|
||||||
tifffolder=os.path.join(outfolder,"tifffolder")
|
tifffolder=os.path.join(outfolder,"舰船","tifffolder")
|
||||||
|
|
||||||
|
|
||||||
h,w=im_data.shape
|
h,w=im_data.shape
|
||||||
|
|
@ -340,9 +331,9 @@ def slicePortDataset(rootname,im_data,src_im_data, im_Geotrans, im_proj, outfold
|
||||||
else:
|
else:
|
||||||
pass
|
pass
|
||||||
# 切片主流程
|
# 切片主流程
|
||||||
binfolder=os.path.join(outfolder,"unit8binfolder")
|
binfolder=os.path.join(outfolder,"港口","unit8binfolder")
|
||||||
pngfolder=os.path.join(outfolder,"pngfolder")
|
pngfolder=os.path.join(outfolder,"港口","pngfolder")
|
||||||
tifffolder=os.path.join(outfolder,"tifffolder")
|
tifffolder=os.path.join(outfolder,"港口","tifffolder")
|
||||||
h,w=im_data.shape
|
h,w=im_data.shape
|
||||||
nextH,ht=getNextSliceNumber(h,portsliceSize,BlockOverLayer)
|
nextH,ht=getNextSliceNumber(h,portsliceSize,BlockOverLayer)
|
||||||
nextW,wt=getNextSliceNumber(w,portsliceSize,BlockOverLayer)
|
nextW,wt=getNextSliceNumber(w,portsliceSize,BlockOverLayer)
|
||||||
|
|
@ -371,18 +362,25 @@ def slicePortDataset(rootname,im_data,src_im_data, im_Geotrans, im_proj, outfold
|
||||||
|
|
||||||
|
|
||||||
def stretchSliceProcess(infilepath, outfolder,portfilestr, strechmethod):
|
def stretchSliceProcess(infilepath, outfolder,portfilestr, strechmethod):
|
||||||
binfolder=os.path.join(outfolder,"unit8binfolder")
|
shipbinfolder=os.path.join(outfolder,"舰船","unit8binfolder")
|
||||||
pngfolder=os.path.join(outfolder,"pngfolder")
|
shippngfolder=os.path.join(outfolder,"舰船","pngfolder")
|
||||||
tifffolder=os.path.join(outfolder,"tifffolder")
|
shiptifffolder=os.path.join(outfolder,"舰船","tifffolder")
|
||||||
allpngfolder = os.path.join(outfolder, "allpngfolder")
|
allpngfolder = os.path.join(outfolder, "allpngfolder")
|
||||||
if not os.path.exists(binfolder):
|
|
||||||
os.makedirs(binfolder)
|
portbinfolder=os.path.join(outfolder,"港口","unit8binfolder")
|
||||||
if not os.path.exists(pngfolder):
|
portpngfolder=os.path.join(outfolder,"港口","pngfolder")
|
||||||
os.makedirs(pngfolder)
|
porttifffolder=os.path.join(outfolder,"港口","tifffolder")
|
||||||
if not os.path.exists(tifffolder):
|
|
||||||
os.makedirs(tifffolder)
|
existOrCreate(shipbinfolder)
|
||||||
if not os.path.exists(allpngfolder):
|
existOrCreate(shippngfolder)
|
||||||
os.makedirs(allpngfolder)
|
existOrCreate(shiptifffolder)
|
||||||
|
existOrCreate(allpngfolder)
|
||||||
|
existOrCreate(portbinfolder)
|
||||||
|
existOrCreate(portpngfolder)
|
||||||
|
existOrCreate(porttifffolder)
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
im_proj, im_Geotrans, im_data=read_tif(infilepath)
|
im_proj, im_Geotrans, im_data=read_tif(infilepath)
|
||||||
src_im_data=im_data*1.0
|
src_im_data=im_data*1.0
|
||||||
im_data = DataStrech(im_data,strechmethod) # 拉伸
|
im_data = DataStrech(im_data,strechmethod) # 拉伸
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue