记2014年数学建模国赛

下面给出，嫦娥三号在粗、精避过程中的算法，python 2.7

#coding=gbk
'''
Created on 2014-9-10

@author: Administrator
'''
from libtiff import TIFF
import numpy
import math
#filename='C:\3.tif'                        #导入附件三文件
filename='C:\4.tif'                          #导入附件四文件
tif = TIFF.open(filename, mode='r')
image = tif.read_image()


def seq(p,a):                               #计算每个区域的标准差函数
    sum = 0.0
    sum1 =0.0
    
    narray=numpy.array(p)
    sum = narray.sum()
    ave = sum/10000.0           
    
    for i in range(100):
        for j in range(100):
            sum1 = sum1 + int(p[i][j])*int(p[i][j])
    s = (sum1-10000.0*(ave**2))/10000.0
    if math.sqrt(s)<5:                    #人为设定的筛选精度
        print '第%d个块的标准差为:%f' %(a,math.sqrt(s))

#size=23                                #2.4km处大图分割参数
size=10                                 #100m处大图分割参数
a = 0
for i in range(size):                           #迭代求解
    for j in range(size):
        q = image[i*100:(i+1)*100,j*100:(j+1)*100]
        a = a+1
        seq(q,a)
        
tif.close()