00001
00002
00003
00004
00005
00006
00007
00008
00009
00010 static const char docstring[] =
00011 "hmi_segment driver for HMI segmentation\n"
00012 "\n"
00013 " [y,s,post,nclean]=hmi_segment(xm,xp,edge,iter,T,beta,alpha,geom,rho,m1,...)\n"
00014 " * Integrated routine for deriving HMI segmentations. Uses\n"
00015 " models (m1,m2,...), plus a magnetogram-photogram pair (xm,xp),\n"
00016 " to deduce an integer labeling. Besides images and models,\n"
00017 " it also requires some disk parameters and labeling smoothness\n"
00018 " parameters.\n"
00019 " * The defaults for iter and T are inherited from mrf_segment_wts,\n"
00020 " and should be looked up there.\n"
00021 " * The posterior probability `post' of the mrf fit can be returned;\n"
00022 " this can quickly verify repeatability.\n"
00023 " * The extreme limb can be a problem for labeling; `edge' helps with\n"
00024 " this. We set the last edge(1) pixels to the value edge(2), and\n"
00025 " return the number of such pixels modified, optionally, in nclean.\n"
00026 " Typically edge=[3 1]; both can be real-valued.\n"
00027 "\n"
00028 " Inputs:\n"
00029 " real xm(m,n)\n"
00030 " real xp(m,n)\n"
00031 " real edge(2)\n"
00032 " int iter(1) or (2)\n"
00033 " real T[0] or [1] or [2] or [3] or [4]\n"
00034 " real beta[1] or [K,K]\n"
00035 " real alpha[K] or [0] = []\n"
00036 " real geom(5)\n"
00037 " real rho\n"
00038 " real m1(l,k1)\n"
00039 " ...\n"
00040 " real mR(l,kR)\n"
00041 "\n"
00042 " Outputs:\n"
00043 " int y(m,n)\n"
00044 " opt real s(R,nS)\n"
00045 " opt real post\n"
00046 " opt int nclean\n"
00047 "\n"
00048 " See Also: makemrfdiscwts, mrf_segment_wts, mixNprobNd,\n"
00049 " clean_edge_label, roi_stats_mag\n"
00050 "\n"
00051 " turmon oct 2009, june 2010, june 2011\n"
00052 "\n"
00053 "";
00054
00055
00056