《西安邮电大学图像处理技术实验报告.docx》由会员分享，可在线阅读，更多相关《西安邮电大学图像处理技术实验报告.docx（74页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、西安邮电大学图像处理技术实验报告西安邮电高校图像处理技术试验报告 本文关键词：西安，图像处理，邮电高校，试验，报告西安邮电高校图像处理技术试验报告 本文简介：装订线西安邮电高校通信与信息工程学院图像处理技术课内试验报告（2022/2022学年第1学期）学生姓名：专业班级：学号：指导老师：西安邮电高校通信与信息工程学院图像处理技术课西安邮电高校图像处理技术试验报告 本文内容：装订线西安邮电高校通信与信息工程学院图像处理技术课内试验报告（2022/2022学年第1学期）学生姓名：专业班级：学号：指导老师：西安邮电高校通信与信息工程学院图像处理技术课内试验过程考核表学生姓名班级/学号担当任务试验室（
2、单位）信工试验室(3#444、436)所在部门信息工程系实施时间2022年8月29日2022年12月26日详细内容试验1：Matlab图像处理基础试验试验2：图像空间坐标变换试验试验3：直方图均衡化图像增加试验试验4：空域滤波图像增加试验试验5：彩色图像增加试验试验6：图像变换及频谱分析试验试验7：图像低通与高通滤波试验试验8：交互消噪及光照校正试验指导老师（师傅）姓名李大湘职务或职称副教授指导老师（师傅）对学生的评价学习看法仔细一般不仔细学习纪律全勤间或缺勤常常缺勤实践实力很强一般较差指导老师（师傅）对学生专业学问或编程实践实力等状况的评语看法指导老师（师傅）签字*年*月*日西安邮电高校通信
3、与信息工程学院图像处理技术课内试验成果鉴定表学生姓名专业班级/学号进行时间2022年8月29日2022年12月26日考核内容标准平常成果（30分）课堂出勤状况、课堂纪律及学习看法（25分）遵守试验室的规章制度状况（5分）试验实力（40分）分析、解决实际问题实力（10分）文件检索或软硬件设计实力（10分）试验完成状况（20分）实习报告（30分）语言及文字表达力（10分）试验报告结构（10分）书写格式规范（10分）指导老师李大湘实习成果成果对应等级试验1：Matlab图像处理基础试验1一、试验目的1二、试验原理1三、试验内容（调试好的程序，试验结果及分析）2试验2：图像空间坐标变换试验5一、试验目
4、的5二、试验原理5三、试验内容（调试好的程序，试验结果及分析）6试验3：直方图均衡化图像增加试验9一、试验目的9二、试验原理9三、试验内容（调试好的程序，试验结果及分析）10试验4：空域滤波图像增加试验14一、试验目的14二、试验原理14三、试验内容（调试好的程序，试验结果及分析）14试验5：彩色图像增加试验19一、试验目的19二、试验原理19三、试验内容（调试好的程序，试验结果及分析）20试验6：图像变换及频谱分析试验25一、试验目的25二、试验原理25三、试验内容（调试好的程序，试验结果及分析）26试验7：图像低通与高通滤波试验.31一、试验目的31二、试验原理31三、试验内容（调试好的程
5、序，试验结果及分析）31试验8：交互消噪及光照校正试验38一、试验目的38二、试验原理38三、试验内容（调试好的程序，试验结果及分析）38试验总结及心得体会42指导老师评语：4344试验1：Matlab图像处理基础试验一、试验目的1.驾驭Matlab图像文件的打开与显示等基本方法；2.驾驭图像数据类型转换及图像类型转换；3.熟识图像矩阵的基本操作。二、试验原理1.MATLAB中图像数据的读写及显示：imread、imshow(1)imread：imread函数用于读入各种图像文件，其一般的用法为语法：X，MAP=imread(filename,fmt)其中，X，MAP分别为读出的图像数据和颜色
6、表数据，fmt为图像的格式，filename为读取的图像文件（可以加上文件的路径）。(2）imwrite:imwrite函数用于输出图像，其语法格式为：语法：imwrite(X,map,filename,fmt)其中，fmt指定的格式将图像数据矩阵X和调色板map写入文件filename。(3)imfinfo：imfinfo函数用于读取图像文件的有关信息，其语法格式为语法：imfinfo(filename,fmt)imfinfo函数返回一个结构info，它反映了该图像的各方面信息，其主要数据包括：文件名（路径）、文件格式、文件格式版本号、文件的修改时间、文件的大小、文件的长度、文件的宽度、每个
7、像素的位数、图像的类型等。(4)imshow：imshow函数是最常用的显示各种图像的函数，语法：imshow(X,map)其中X是图像数据矩阵,map是其对应的颜色矩阵，若进行图像处理后不知道图像数据的值域可以用代替map。须要显示多幅图像时，可以运用figure语句，它的功能就是重新打开一个图像显示窗口。2.程序流程限制语句：for循环结构、while循环结构、if-else-end分支结构(1)for循环结构语法：fori初值：增量：终值语句1语句nend说明：i初值：终值，则增量为1。初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑(2)while循环结构语法：wh
8、ile逻辑表达式循环体语句end(3)if-else-end分支结构if表达式1语句1elseif表达式2(可选)语句2else(可选)语句3endend说明：i.if结构是一个条件分支语句，若满意表达式的条件，则往下执行；若不满意则跳出if结构。iielseif表达式2与else为可选项，这两条语句可依据详细状况取舍。iii留意：每一个if都对应一个end，即有几个if，记就应有几个end。三、试验内容（调试好的程序，试验结果及分析）1、图像的打开与显示。在matlab吩咐窗口执行：程序;f=imread(e:/imdata/T1.jpg);imshow(f)；图1-12、图像的信息查询。在
9、matlab吩咐窗口执行：info=imfinfo(ngc6543a.jpg)图1-23、图像的存储。在matlab吩咐窗口执行：图1-34、彩色图像灰度化编程试验。在matlab打开文本编辑器，编写Lab1_Test1.m文件：程序代码：clearall;closeall;clcf=imread(e:/imdata/onion.png);figure,subplot(1,2,1),imshow(f),title(原图);g=rgb2gray(f);subplot(1,2,2),imshow(g),title(灰度图像);图1-45、图像数值算术运算编程试验。在matlab打开文本编辑器，编写
10、Lab1_Test2.m文件，实现：（1）将每个像素加上50；（2）255减每个像素的值：6、图像逻辑运算编程试验。在matlab用如下语句构造二幅“二值图像”，自己打开文本编辑器，编写Lab1_Test3.m文件，实现：（1）统计图像A与B中非0像素的数量；（2）将图像A与B“相或”“相与”，得到新的一幅图像C；视察比较运算结果。试验2：图像空间坐标变换试验一、试验目的1.驾驭Matlab图像文件的打开与显示等基本方法；2.驾驭图像数据类型转换及图像类型转换；3.熟识图像矩阵的基本操作。二、试验原理图像的几何变换由两个算法来实现：空间变换算法和图像插值算法。空间变换算法采纳数学函数形式描述输
11、入、输出图像对应像素间的空间关系。设原图像的(x,y)点变换后移动到新图像的(x1,y1)点，(x,y)和(x1,y1)的空间变换关系可由下式描述：其中，T为3*3的变换矩阵。1、图像平移：初始坐标为的点经过平移，坐标变为，两点之间的关系为：，以矩阵形式表示为：2、图像的镜像变换：是以图象垂直中轴线或水平中轴线交换图像的变换，分为垂直镜像变换和水平镜像变换，两者的矩阵形式分别为：3、图像缩小和放大（缩放）：变换矩阵相同：当时，图像缩小；时，图像放大。4、图像旋转：定义为以图像中某一点为原点以逆时针或顺时针方向旋转肯定角度。其变换矩阵为：该变换矩阵是绕坐标轴原点进行的，假如是绕一个指定点（a,b
12、）旋转，则现要将坐标系平移到该点，进行旋转，然后再平移回到新的坐标原点。3、试验内容（调试好的程序，试验结果及分析）1、试用Matab自带的imresize（）函数，实现图像缩放，imresize函数常用的语法格式如下：(1)用不同的插值方法实现图像缩放Im=imresize(I,m,method);参数I：待缩放的图像，由imread函数读入。参数m：图像放大倍数，m1)Img=rgb2gray(Img);endfigure;subplot(2,2,1);imshow(Img);title(原图像);Hist=zeros(1,256);fori=1:heightforj=1:widthHis
13、t(Img(i,j)+1)=Hist(Img(i,j)+1)+1;endendsubplot(2,2,2);bar(Hist);title(原图像的直方图);p=zeros(1,256);fori=1:256p(i)=Hist(i)/(height*width*1.0);endc=zeros(1,256);c(1)=p(1);fori=2:256c(i)=c(i-1)+p(i);endc=uint8(255.*c+0.5);fori=1:heightforj=1:widthI(i,j)=c(Img(i,j)+1);endendsubplot(2,2,3);imshow(I);title(均衡化
14、的图像);subplot(2,2,4);imhist(I);title(新图像的直方图);图3-3试验4：空域滤波图像增加试验一、试验目的1.驾驭Matlab图像文件的打开与显示等基本方法；2.驾驭图像数据类型转换及图像类型转换；3.熟识图像矩阵的基本操作。二、试验原理1、空间滤波机理：就是在待处理图像中逐点地移动掩模。在每一点(x,y)处，滤波器在该点的响应通过事先定义的关系来计算。对于线性空间滤波，其响应由滤波器系数与滤波掩模扫过区域的相应像素值的乘积之和给出。2、均值滤波器：运用滤波器模板确定的平均灰度值代替图像中的每个像素的灰度值，这种处理的结果降低了图像灰度的“尖锐”改变。由于典型的
15、的随机噪声有灰度级的急剧改变组成，因此，常见的均值滤波器处理应用就是降低噪声。但是由于图像边缘也是由图像灰度尖锐改变带来的特性，所以均值滤波处理还是存在着不希望有的边缘模糊的负面效应。3、加权平均滤波器：所谓加权是指用不同的系数乘以像素，使得一些像素的重要性（权重）比另一些像素的重要性更大。假如给予中心点最高权重，然后随着距中心点距离的增加而减小系数的加权权重，其目的是在平滑处理中试图降低模糊。也可以选择其他权重来达到所需的目的。4、中值滤波器：所谓中值滤波就是图像中每个像素为中心，把它一个邻域中全部像素取出来，由大到小进行排序，中间位置那个值为中值，用它代替图像中的每个像素的灰度值。5、锐化
16、空间滤波器：运用二阶微分（拉普拉斯算子）进行图像锐化，由于拉普拉斯算子是一种微分算子，因此其应用强调的是图像的突变，并不强调灰度级缓慢改变的区域。将原图像和拉普拉斯图像叠加在一起，可以复原背景特性并保持拉普拉斯锐化处理的效果，假如运用的定义是具有负的中心系数的拉普拉斯滤波模三、试验内容（调试好的程序，试验结果及分析）1、均值滤波试验。在Matlab中，建立M文件Lab4_Test1.m，编写如下代码，利用Matlab自带的函数对噪声图像进行滤波处理，且分析对比试验结果。备注：方法一：用filter2函数进行滤波处理。方法二：用conv2卷积函数进行滤波处理。程序4-1（1）functionLa
17、b4_Test1_Main()clc;closeall;clearall;I=imread(E:/imdata/lena.jpg);I=rgb2gray(I);figure;subplot(1,3,1);imshow(I);title(原图像);M=imnoise(I,gaussian,0,0.05);subplot(1,3,2);imshow(M);title(加噪图像);J=filter2(fspecial(average,3),M)/255;subplot(1,3,3);imshow(J);title(滤波后图像);程序4-1（2）clc;closeall;clearall;I=imre
18、ad(E:/imdata/lena.jpg);I=rgb2gray(I);M=imnoise(I,gaussian,0,0.05);H1=ones(3)/9;H2=ones(7)/49;M=double(M);G1=conv2(M,H1,same);G2=conv2(M,H2,same);figure;subplot(2,2,1);imshow(I);title(原图像);subplot(2,2,2);imshow(M,);title(添加噪声的图像);subplot(2,2,3);imshow(G1,);title(3*3均值滤波的图像);subplot(2,2,4);imshow(G2,)
19、;title(7*7均值滤波的图像);图4-12、自编程序均值滤波试验。在Matlab中，建立M文件Lab4_Test2.m，编写如下代码，对噪声图像进行均值滤波处理，且分析对比试验结果。functionLab4_Test2_Main()clc;closeall;clearall;Img=imread(E:/imdata/lena.jpg);%打开图像K=5;%邻域半径，自己改此值，对比滤波效果II=MyMeanFilter(Img,K);%调用自编直方均衡化disp(RunOK)functionI=MyMeanFilter(Img,K)%直方图均衡化height,width,dim=size
20、(Img);if(dim1)Img=rgb2gray(Img);%若打开的为彩色图像，则先转化成灰度图像endfiguresubplot(1,3,1),imshow(Img),title(原图像);%显示原始图像M=imnoise(Img,gaussian,0,0.05);subplot(1,3,2),imshow(M),title(加噪图像);I=M;R=floor(K/2);%邻域宽度的一半forj=R+1:height-Rfori=R+1:width-Rblock=Img(j-R:j+R,i-R:i+R);t=sum(block(:);%求和t=t/(K*K);%平均值I(j,i)=t;
21、endendI=uint8(I);%数据类型由double变成uint8subplot(1,3,3),imshow(I),title(加噪图像);图4-2（1）图4-2（2）3、中值滤波试验。在Matlab中，建立M文件Lab4_Test3.m，编写如下代码，实现中值滤波功能对噪声图像进行滤波处理，且分析对比试验结果。functionLab4_Test3_Main()clc;closeall;clearall;Img=imread(e:/T1.jpg);%打开图像K=5;%邻域半径，自己改此值，对比滤波效果II=MyMedfilt2(Img,K);%调用自编直方均衡化disp(RunOK)fu
22、nctionI=MyMedfilt2(Img,K)%直方图均衡化height,width,dim=size(Img);if(dim1)Img=rgb2gray(Img);%若打开的为彩色图像，则先转化成灰度图像endfiguresubplot(1,3,1),imshow(Img),title(原图像);%显示原始图像M=imnoise(Img,saltsubplot(1,3,2),imshow(M),title(加噪图像);I=M;R=floor(K/2);%邻域宽度的一半Pos=floor(K*K/2)+1;%中值位置forj=R+1:height-Rfori=R+1:width-Rbloc
23、k=Img(j-R:j+R,i-R:i+R);sequece=sort(block(:);%排序media=sequece(Pos);I(j,i)=media;endendI=uint8(I);%数据类型由double变成uint8subplot(1,3,3),imshow(I),title(加噪图像);图4-34、锐化滤波试验。以下图所示模板为例，在Matlab中，建立M文件Lab4_Test4.m，参照上述编程风格，自己编写程序实现“拉普拉斯图像锐化”功能，且分析对比试验结果。0101-41010functionchap3_Laplacesharp()clearall;closeall;c
24、lc;I=imread(e:/T1.jpg);I=rgb2gray(I);figure;subplot(1,2,1);imshow(I);title(原图);H,W=size(I);I=double(I);G=I;fory=2:H-1forx=2:W-1d2=I(y-1,x)+I(y+1,x)+I(y,x-1)+I(y,x+1)-4*I(y,x);t=I(y,x)+abs(d2);if(t255)t=255;endG(y,x)=t;endendsubplot(1,2,2);imshow(uint8(G);title(锐化后结果);disp(程序结束);图4-4试验5：彩色图像增加试验一、试验目
25、的1.驾驭Matlab图像文件的打开与显示等基本方法；2.驾驭图像数据类型转换及图像类型转换；3.熟识图像矩阵的基本操作。二、试验原理1、伪彩色图像增加：伪彩色(pseudo-color)图像的每个像素值事实上是一个索引值或代码，该代码值作为色调查找表CLUT(ColorLook-UpTable)中某一项的入口地址，依据该地址可查找出包含实际R、G、B的强度值。这种用查找映射的方法产生的色调称为伪彩色。用这种方式产生的色调本身是真的，它不肯定反映原图的色调。伪彩色增加的方法主要有：密度分割、彩色变换和频率域伪彩色增加三种。1)密度分割法密度分割法又称为灰度分割法，是伪彩色处理技术中最基本，最简
26、洁的方法。设一幅灰度图像f(x,y),可以看成是坐标(x,y)一个密度函数。把此图像的灰度分成若干等级，即相当于用一些和坐标平面平行的平面切割此密度函数。通常的作法是：用一个平行于图像坐标平面XY的平面去切割图像亮度函数；对每一个输入灰度值，假如它在切割灰度值之上就给予某一种颜色，假如它在之下就给予另一种颜色。2)灰度变换法对每个原始图中像素的灰度值用三个独立变换来处理，从而将不同的灰度映射为不同的彩色变换后原始图中灰度值偏小的像素将主要呈现绿色，而灰度值偏大的像素主要呈现红色。2、真彩色图像增加：真彩色是指图像中的每个像素值都分成R、G、B三个基色重量，每个基色重量干脆确定其基色的强度，这样
27、产生的色调称为真彩色。也有称自然物体的彩色叫真彩色。真彩色图像就是能真实反映自然物体原来颜色的图像。处理策略1）将一幅彩色图像看作三幅重量图像的组合体，在处理过程中先对每幅图像（根据对灰度图像处理的方法）单独处理，再将处理结果合成为彩色图像2）将一幅彩色图像中的每个像素看作具有三个属性值，即像素属性现在为一个矢量，需利用对矢量的表达方法进行处理，可进行模板处理。三、试验内容（调试好的程序，试验结果及分析）1、彩色图像通道分别与颜色空间转化试验。在Matlab中，建立M文件Lab5_Test1.m，编写如下代码，打开RGB彩色图像，将转化到HSV空间，且分别显示每个通道。functionLab5
28、_Test1()clc;closeall;clearall;%说明：打开RGB图像，分别三个通道分别显示出来Img=imread(e:/imdata/lena.jpg);%打开图像figure,subplot(1,4,1),imshow(Img),title(RGB原图像);%分别三个通道R=Img(:,:,1);%分别取三维数组的一维，得到红绿蓝三个重量G=Img(:,:,2);%为RGB。B=Img(:,:,3);subplot(1,4,2),imshow(R),title(R通道图像);subplot(1,4,3),imshow(G),title(G通道图像);subplot(1,4,4
29、),imshow(B),title(B通道图像);%说明：RGB空间转化到HSV空间hsv=rgb2hsv(Img);%分别三个通道H=hsv(:,:,1);%分别取三维数组的一维，得到红绿蓝三个重量S=hsv(:,:,2);%为RGB。V=hsv(:,:,3);figure,subplot(1,4,1),imshow(hsv),title(HSV原图像);subplot(1,4,2),imshow(H),title(H通道图像);subplot(1,4,3),imshow(S),title(S通道图像);subplot(1,4,4),imshow(V),title(V通道图像);disp(程
30、序胜利结束)图5-1(1)图5-1（1）2.伪彩色图像增加试验亮度切分法在Matlab中，建立M文件Lab5_Test2.m，编写如下代码，采纳亮度切分法将一幅灰度图像转化成伪彩色图像。3.伪彩色图像增加试验灰度变换法在Matlab中，建立M文件Lab5_Test3.m，编写如下代码，采纳三种不同的灰度变换就法，将一幅灰度图像转化成伪彩色图像。functionLab5_Test2()clearall;closeall;clc%I=imread(1.jpg);I=imread(e:/imdata/lena.bmp);figure;subplot(1,2,1),imshow(I),title(a)
31、原始图像);M,N=size(I);Img=double(I);%把数据变成双精度型Img=Img./255;%数据范围归到01区间%采纳三种不同的变换对数据进行变换forx=1:Mfory=1:N%第一种变换R(x,y)=255*log2(1+Img(x,y);%其次种变换t=2*Img(x,y);if(t1)Img=rgb2gray(Img);%若打开的为彩色图像，则先转化成灰度图像endfigure,subplot(1,3,1),imshow(Img),title(原图像);%2图像FT变换及频谱F=fft2(Img);%计算二维傅里叶变换F=fftshift(F);%将原零点移到频谱图
32、的中心%显示变换后的频谱图subplot(1,3,2),imshow(log(abs(F)+1),010),title(频谱);%3图像IFT变换F1=ifftshift(F);%系数逆中心化f1=ifft2(F1);%逆傅里叶变换复原图像subplot(1,3,3),imshow(uint8(f1),title(IFFT复原的图像);%主程序到此结束,图6-12、图像的DCT变换试验。在Matlab中，建立M文件Lab6_Test2.m，编写如下代码，打开灰度图像，进行DCT变换，显示频谱并进行IDCT变换复原图像。functionLab6_Test2()clc;closeall;clear
33、all;%说明：打开进行DCT变换，显示频谱，再进行IDCT变换复原图像%1打开图像Img=imread(E:/imdata/lena.jpg);%打开图像H,W,dim=size(Img);if(dim1)Img=rgb2gray(Img);%若打开的为彩色图像，则先转化成灰度图像endfigure,subplot(1,3,1),imshow(Img),title(原图像);%2图像DCT变换及频谱F=dct2(Img);%计算二维傅里叶变换%显示变换后的频谱图subplot(1,3,2),imshow(log(abs(F)+1),010),title(DCT频谱);%3图像IDCT变换f1
34、=idct2(F);subplot(1,3,3),imshow(uint8(f1),title(IDCT复原的图像);图6-23、图像FT旋转定理验证明验。在Matlab中，建立M文件Lab6_Te篇2：西安市河道采砂管理方法西安市河道采砂管理方法 本文关键词：西安市，河道，管理方法，采砂西安市河道采砂管理方法 本文简介：西安市人民政府令第56号西安市河道采砂管理方法已经2022年4月20日市人民政府第78次常务会议通过，现予公布，自2022年8月1日起施行。市长孙清云二00五年六月十一日西安市河道采砂管理方法第一条为加强河道的整治和管理，维护河势和河床稳定，保障防洪工程、铁路、马路、桥梁等设
35、施的平安和行洪畅通，西安市河道采砂管理方法 本文内容：西安市人民政府令第56号西安市河道采砂管理方法已经2022年4月20日市人民政府第78次常务会议通过，现予公布，自2022年8月1日起施行。市长孙清云二00五年六月十一日西安市河道采砂管理方法第一条为加强河道的整治和管理，维护河势和河床稳定，保障防洪工程、铁路、马路、桥梁等设施的平安和行洪畅通，合理开发利用河道砂石资源，依据中华人民共和国水法、陕西省河道管理条例等有关法律、法规，结合本市实际，制定本方法。其次条本市行政区域内的河道采砂管理适用本方法。第三条本方法所称河道采砂，是指在河道管理范围内采运砂石、取土和淘金（包括淘取其他金属和非金属
36、）等活动。第四条市水行政管理部门是本市河道采砂的主管机关。区、县水行政管理部门根据河道管理权限，负责本行政区域内河道采砂的管理和监督检查工作，并接受市水行政管理部门的监督与指导。第五条河道采砂实行统一管理与分级管理相结合以及依法开采、合理利用和统筹兼顾、实行总量限制的原则，根据有利于发展经济和维系优良水生态环境的要求，保障河势稳定、行洪通畅及堤防、铁路、马路、桥梁、管道、输电线路、通信电缆、水文观测设施及临河工程设施的平安。第六条水行政管理部门应依据所管辖河道、河段的河势现状和砂、石、土料状况，结合河道防洪规划、整治规划，对河道采砂进行统一规划，划定可采区、禁采区，设定禁采期，限定开采深度和年
37、度开采总量。对边界河道采砂，应根据河道中心线和民政部门划定的行政边界实行属地管理。第七条本市渭河段的河道采砂规划由市水行政管理部门组织编制，报上级水行政管理部门审批；其他河道的采砂规划由河道管辖地水行政管理部门编制，报市水行政管理部门审批。河道采砂规划一经批准，应当严格执行，确需修改的，应当经原批准机关同意。河道采砂规划涉及铁路、马路、电力、通信、管道等设施平安爱护范围的，编制时应征求相关管理部门的看法。第八条水行政管理部门可以依据河道或河段砂源补给状况、河道采砂对防洪工程及相关建筑设施的影响，对采砂已严峻超限，河床严峻下切，河道两岸地下水位严峻下降的河道、河段，提出全面禁止采砂的看法和方案，
38、经同级人民政府批准后发布公告，并报上一级水行政管理部门备案。第九条河道采砂实行许可证制度。从事河道采砂的单位和个人应当持有河道采砂许可证。对申请及审批发放河道采砂许可证的基本条件和基本程序，水行政管理部门应当进行公示。从事河道采砂的单位和个人应当向有管辖权的水行政管理部门提出河道采砂申请。水行政管理部门自接到河道采砂单位或个人申请之日起20个工作日内应完成审批程序，符合规划要求和开采条件的，发给河道采砂许可证；不符合条件的，应予书面答复并说明理由。对符合规划要求和开采条件的河段，水行政管理部门可实行公开招标的方式确定采砂人，也可公开拍卖采砂权。第十条河道采砂许可证每年汛期结束后发放，有效期限不
39、超过一年。不得伪造、转让、涂改、买卖、出租、出借河道采砂许可证。第十一条从事经营性开采的，应根据有关规定到当地工商、税务等相关管理部门办理有关手续。第十二条经批准在河道管理范围内采砂的单位或个人，必需遵守下列规定：（一）每月月底前向批准开采的水行政管理部门或其托付的河道管理单位报送本月采砂状况及下月开采安排；（二）严格根据批准的开采地点、开采范围、开采深度、作业方式、开采时间、砂石料堆放地点作业；（三）随时清除或复平弃料堆体；（四）不得损坏测量标记、水文观测设施、照明报警设施、通信电缆、宣扬牌、界桩、里程桩、堤顶路面、护堤护林设施和河道防护林等；（五）不得阻碍河道行洪和防洪抢险。第十三条对河道
40、管理范围内的防洪工程、水利枢纽、水文观测设施、护堤地、护岸地以及铁路、马路、桥梁、管道、输电线路、通信电缆等工程设施，各级水行政管理部门应会同相关管理部门依法划定平安爱护范围，并设立明确标记。禁止在防洪工程、水工程、铁路、马路、桥梁、水库大坝、涵闸、水文观测设施、通信电缆、输气管道、输电线路等工程设施平安爱护范围内采砂。禁止在防洪工程上堆放砂石料。第十四条水行政管理部门或其托付的河道管理单位应定期对所管辖的河道采砂进行检查，对可能影响防洪工程和其他工程设施平安的，水行政管理部门可以变更或限制开采范围、开采量和作业方式，必要时应当责令其停止河道采砂，并收回河道采砂许可证。第十五条汛期内禁止运用采
41、砂机械在河道采砂，放置在河道内的采砂机械和砂石堆体必需在汛期到来前撤出河道管理范围。在汛情紧急状况下，水行政管理部门和防汛指挥机构有权停止河道采砂作业，并清除有碍行洪的砂石料堆体和采砂设备、设施。第十六条经批准在河道采砂的单位或个人应当为其运输车辆修筑越堤路。修筑越堤路的地址和方案，应当经河道所在地水行政管理部门或其托付的河道管理单位审查批准。运砂车辆应根据指定的路途行驶，不得损坏防洪工程、影响防洪抢险。第十七条经批准在河道采砂的单位或个人，应依法向发放河道采砂许可证的水行政管理部门或其托付的河道管理单位缴纳河道砂石资源费。河道砂石资源费征收标准按省物价部门核定的标准执行，不得随意减免、提高。
42、河道砂石资源费全额上缴财政，专项用于河道治理和管理。第十八条河道沿岸村民个人运用非机械设施自采自用少量（10立方米以下）河道砂、石、土料的，可干脆向所在地水行政管理部门或其托付的河道管理单位提出申请，经批准后在指定地点采运，不需办理河道采砂许可证，免征河道砂石资源费。第十九条违反本方法规定，未办理河道采砂许可证擅自由河道管理范围内采砂的，由水行政管理部门责令停止违法行为、暂扣非法采砂机具、没收违法所得，并处10000元以上20000元以下罚款。在禁采区、禁采期采砂的，处20000元以上30000元以下罚款。其次十条违反本方法第十条规定，伪造、转让、涂改、买卖、出租、出借河道采砂许可证的，收缴或
43、吊销其河道采砂许可证，对单位处1000元以下罚款，个人处200元以下罚款；对有违法所得的，没收违法所得，并处10000元以下罚款。其次十一条违反本方法第十二条第（二）、（三）项规定的，由水行政管理部门责令限期改正，并处10000元以上20000元以下罚款。其次十二条违反本方法第十六条规定的，由水行政管理部门责令其改正，情节严峻的，吊销河道采砂许可证，并可处以单位1000元、个人200元罚款。其次十三条违反本方法第十七条规定，不按时缴纳砂石资源费的，每逾期一天加收砂石资源费总额3的滞纳金。其次十四条当事人对行政惩罚确定不服的，可依法申请行政复议或提起行政诉讼。其次十五条河道采砂管理人员玩忽职守、滥用职权、徇私舞弊或因管理不善造成河道采砂秩序混乱和责任事故的，由所在单位或上级主管机关赐予行政处分；造成重大损失或构成犯罪的，由司法机关依