超人琪的博客

原创 非线性函数模型（多项式拟合）预测土壤侵蚀模数

某地区各地理单元的土壤侵蚀模数（x）与土壤含氮量1（y）的数据见（实验1数据中的“题目3”sheet）。画出二者之间的散点图并确定是什么样的相关形式； 试用一个非线性函数模型拟合该地区土壤侵蚀模数（x）与土壤含氮量（y）的关系； 检验该模型的显著性，并预测当土壤侵蚀模数x=8000t/(km2*a)时的土壤含氮量y。1. # coding=utf-8 2. import matplotlib.pyplot as plt 3. from matplotlib.font_manager.

原创 计算偏相关系数和复相关系数

根据4个要素的48个样本数据，计算得到的简单相关系数如下，试计算各级偏相关系数和复相关系数，并对其显著性进行检验。R=r11r12r13r14r21r22r23r24r31r32r33r34r41r42r43r44=10.99540.99040.99890.995410.99260.99230.99940.992610.99770.99890.99230.997711. import math 2. import numpy as np 3. import pandas as pd..

原创 计算相关系数进行显著性检验

某地区粮食产量与受灾面积的历年数据见（实验1数据中的“题目1”sheet），试计算二者之间的相关系数，并对相关系数进行显著性检验（α =0.05）1. import numpy as np 2. 3. x_simple=np.array([251,801,200,409,415,502,314,1101,980,1124]) 4. y_simple=np.array([52,101,65,88,90,98,120,150,140,120]) 5. 6. #计算相关系...

原创 python数据分析——文件操作

文件的打开>>> f1=open('d:\\infile.txt')//读入文件>>> f2=open(r'd:\outfile.txt','w')//写文件、清空文件内容或新建一个文件>>> f3=open('record.dat','wb',0)//a是追加//r+=r+w//w+=w+r/a+=a+r文件相关函数...

原创 AｒｃＧｉｓ提取某一类型的土地

１、图层点击右键打“打开属性表”是灰色的怎么办？并不是所有栅格都是自带属性表的啦工具箱－ｄａｔａ　ｍａｎａｇｅｍｅｎｔ　ｔｏｏｌ－栅格－栅格属性－构建栅格属性表再次右键打开属性表就可以啦２、如何在土地利用类型图里提取某一土地类型呢？我现在是要在一个县里找到冬小麦的种植区域，然后和其他数据进行叠加分析。在矢量图层上右击打开属性表，在属性表左上角单机按属性选择，看看你的土地利用在哪个字段然后选中他＝获取唯一值，然后选你要选的那个就可以了３、如何在ａｒｃｇｉｓ中建立选区这个

原创 python程序调试与异常

# print(10/0)from ch4.TramError import TranErrortry: # print(10/0) a = b + 1 a = 1 + 1 b=3*2 a=int('abc') raise TranError(500001,'tranError')except (NameError,ZeroDivisionError) as e: print(type(e)) print(e)except Exc.

原创 python——面向对象（继承）

import sys, ossys.path.append(os.path.dirname(os.path.abspath(__file__)))from p1 import Personclass Student(Person): def __init__(self,name,height,stno,school): super().__init__(name,height) #super为继承父类的属性 self.stno =stno .

原创 python——面向对象（封装）

class Person: height=180 #类属性，是实例共享的，成员共享 name='王宝强' count=0 def __init__(self,name,height): #height是形参 # print('i am a person object') # return 'hehe' 必须返回的是空值 self.height=height #实例的属性，相当于person.heigh..

原创 python——面向对象

# 面向对象OO# 类和对象# 属性和方法# 继承、封装、# class ClassDemo: #类名开头必须大写# '''类的一些说明'''# pass## '''类的使用——实例化一个类'''# myclass = ClassDemo()# myclass1 = ClassDemo()## print(myclass)#打印的是对象 <__main__.ClassDemo object at 0x000002D1A06B6AF0># print(m.

原创 python函数调用参数

# 函数中的形参和实参def add(a,b): print("a=",a) print("b=",b) result= a+b return resultc=add(3,1)#必要参数A=add(b=3,a=1) #关键字参数c=add(1)#默认值参数print(c)def print_userinfo(name,gender,age,depart): print("我是："+name) print("我是 "+gender+"生").

原创 python——建立函数

#内置函数 自定义函数# print('hello python')# round(1.12389,2)精确到小数点后2位# 内置函数# 1封装性，隐藏代码的细节# 2避免重复代码# 3函数都有一定的功能# 自定义函数# def funcname(paran s):# pass# 1定义一个打印hello字符串的方法# 2定义一个打印指定字符串的方法# 3定义一个计算加法的函数a+b# def Print(s):# print(s)# ms='he.

原创 IDL基础习语法学习一

IDL基础习语法学习（一）;idl中注释符用';';$续航符'$';&同行符'&';IDL>print,'hello,idl' 命令行模式;（编译代码）.compile "D:\Desktop\毕业论文\IDL_practice\firstsample.pro";批处理模式;CD设为当前目录;IDL>cd,'D:\Desktop\毕业论文\IDL_practice';IDL>@firstbatchsample;或者;IDL>@.

原创 结构体总结二

#include <stdio.h>typedef struct{ int num; char name[9]; float score[3];}STU;void show(STU tt){ int i; printf("%d %s : ",tt.num,tt.name); for(i=0; i<3; i++) printf("%5.1f",tt.score[i]); printf("\n");}/*****.

原创 相关系数法影像匹配

基于VS 2015 opencv#include<opencv2/opencv.hpp>#include<iostream>using namespace std;using namespace cv;int main(){ Mat airport, plane; airport = imread("airport.bmp",0); plane = imread("template.bmp", 0); Mat Correlation(airport.ro.

原创 matlab编程实现Forstner算子

clear all;clc;%获取文件名、路径[Filename,Pathname,FilterIndex]=uigetfile({'test1.png','D:\新建文件夹\遥感\数字摄影测量\Forstner算子'},'Select the input image');if FilterIndex pic_orig=imread([Pathname,Filename]);%读取原始标准影像else returnendfigure(1);imshow(pic_orig);.

原创 新建要素和屏幕数字化

新建要素右键Default新建要素——定义要素名称plant屏幕数字化打开校正后的图层——开始编辑——创建要素（例如：面）——沿边界选择面——点击要素右键——编辑属性...

原创 影像配准和定义投影

影像配准：右键点击空白处，导出地理配准工具条地理配准工具中选择添加控制点找到控制点，左键再右键，输入给定的x,y坐标地理配准工具中选择校正，更新显示，全图显示定义投影：右键属性导出数据图片toolbox数据管理工具投影变换——定义投影——选择图片和投影坐标系...

原创 基础应用——专题地图的制作

在布局视图下操作添加标题，图例，指南针，比例尺，格网摸索到小功能，图例可以分组处理更方便

原创 基础应用——属性表的操作

属性表打开新建属性-人口密度属性表连接编辑器——面积/人口=人口密度成图

原创 基础应用——统计噪音范围内的土地利用

原创 结构体总结一

//找出指定编号人员的数据，作为函数值返回，若指定编号不存在，返回编号为空串#include <stdio.h>#include <string.h>#define N 8typedef struct{ char num[10]; int year,month,day ;}STU;/**********found****...

原创 真题五十五、五十六、五十七

struct student { long sno; char name[10]; float score[3];};void fun(struct student a[], int n){/**********found**********/ struct student t;//结构体类型，参考形参 int i, j;/**********found...

原创 真题五十二、五十三、五十四

//将每个人信息作为数据块写入到二进制文件中#include <stdio.h>#include <stdlib.h>#define N 5typedef struct{ int num; char name[10]; char tel[10];}STYPE;void check();/**********f...

原创 字符串前后*的删减问题

//保留字符串前面n个*号void fun( char *a, int n ){ int i = 0; int k = 0; char *p, *t; p = t = a; while (*t == '*') { k++; t++; } if (k > n)//p的前部保留n各星号，后部依次存入 { while (*p) { a[i] = *...

原创 链表题型总结四

//查找数据中值为ch的结点，返回节点在链表中的顺序#include <stdio.h>#include <stdlib.h>#define N 8typedef struct list{ int data; struct list *next;} SLIST;SLIST *creatlist(char *);voi...

原创 真题四十九、五十、五十一

//最大值放第一个数，最小值放第二个数，次小值放第三个数void fun(int a[], int n){ int i, j, max, min, px, pn, t;/**********found**********/ for (i=0; i<n-1; i+=2) { max = min = a[i]; px = pn = i;/*******...

转载 双线性内插实现图像旋转

谢谢这位优秀的博主：https://blog.csdn.net/wuminmatlab/article/details/43305847?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522158411238619724839200410%2522%252C%2522scm%2522%253A%252220140713.1300568...

原创 真题四十六、四十七、四十八

//字符串变数字long fun ( char *p){ int k = 0; int flag; int a[20] = { 0 }; if (*p == '-') flag = -1; else flag = 1; for (int i = 0; i < strlen(p); i++) { if (p[i] >= '0'&&p[i]...

原创 真题四十三、四十四、四十五

//逆序存放字符串，后面加上原来的字符串//法一：void fun (char *s, char *t){/************found************/ int i; int sl = strlen(s)-1; for (i=0; i<=sl; i++)/************found************/ t...

原创 ３０Ｍ　DEM数据处理成图

最近写论文要用到DEM数据，顺便整理了以下。数据来源：NASA SRTM1 v3.0由于数据是ｈｇｔ文件类型，需要进行栅格数据转换步骤为：打开ＡＲＣＧＩＳ点击目录下Ｄａｔａ　ｍａｎａｇｅｍｅｎｔ　ｔｏｏｌｓ——栅格——栅格数据集——镶嵌ｈｇｔ文件放到输入栅格，目标栅格自己定义名称然后进行矢量裁剪点击目录下Ｄａｔａ　ｍａｎａｇｅｍｅｎｔ　ｔｏｏｌｓ——栅格—...

原创 真题四十、四十一、四十二

//数字字符顺序后移，经常考，要牢记于心char *fun(char *s){ int i, j, k, n; char *p, *t; n=strlen(s)+1; t=(char*)malloc(n*sizeof(char)); p=(char*)malloc(n*sizeof(char)); j=0; k=0; for(i=0; i<n; i++)...

原创 真题三十七、三十八、三十九

注意细节修改题里，（*s)++改成s++，一定要把括号去掉while(最好是填字符，单向链表填p，字符串数组填*p)//判断字符串是否结束//将字符串中标为奇数的字符取出，从小到大排列void fun(char *s, char *p){ int i, j, n, x, t; n=0; for(i=0; s[i]!='\0'; i++) n++; f...

原创 python语法学习六——异常

try-except语句>>> try: num1=int(input('1:')) num2=int(input('2:')) print(num1/num2)except ValueError: print('digit:')>>> try: num1=int(input('1:')) num2=int(input('2:')) p...

原创 python语法学习五——函数库

math库>>> import math>>> dir(math)['__doc__', '__loader__', '__name__', '__package__', '__spec__', 'acos', 'acosh', 'asin', 'asinh', 'atan', 'atan2', 'atanh', 'ceil', 'copysi...

原创 python语法学习四——函数

break语句>>> i=1>>> while i%3: print(i,end='') if i>=10: break i+=1 12continue语句自定义函数#输出素数from math import sqrtdef isprime(x): if x == 1: return F...

原创 真题三十四、三十五、三十六

//函数定义；结构体变量成员的表示方法#include <stdio.h>#include <string.h>struct student { long sno; char name[10]; float score[3];};/**********found**********/struct student fun(stru...

原创 链表题型总结三

//将不带头结点的单向链表逆置#include <stdio.h>#include <stdlib.h>#define N 5typedef struct node { int data; struct node *next;} NODE;/**********found**********///重点掌握NODE * fun(N...

原创 真题三十一、三十二、三十三

//在每个数字字符后插入*void fun(char *s){ int i, j, n; for(i=0; s[i]!='\0'; i++)/**********found**********/ if(s[i]>='0' && s[i]<='9') { n=0;/**********found**********/ ...

原创 真题二十八、二十九、三十

//复习switch结构的用法#include <stdio.h>void fun( int *s, int n ){/**********found**********/ int i, one=0, zero=0 ; for(i=0; i<n; i++)/**********found**********/ switch( s[i] )...

原创 链表题型总结二

//将a,b,c各结点的数据域中存入字符#include <stdio.h>typedef struct list{ char data; struct list *next;} Q;void fun( Q *pa, Q *pb, Q *pc){ Q *p;/**********found**********/ pa->next=pb;...