Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1264 1711 (N47W34 L=134)  1505    279        3 (1231/1619-1322/1784)  34.5 W27-W44 N40-N52
      2  930 1499 (N37E00 L=100)   303     50        1 ( 900/1469- 960/1517) -14.0 E02-W02 N35-N38
      3  990 1427 (N31W03 L=103)   264     73        0 ( 964/1396-1016/1458) -68.2 W01-W04 N28-N33
      4  719 1384 (N32E17 L= 83)   401    131        0 ( 691/1332- 748/1439) -60.3 E18-E14 N28-N35
      5 1562 1221 (N10W38 L=138)   432     94        2 (1517/1181-1595/1266) -71.6 W35-W40 N07-N13
      6 1233 1183 (N13W15 L=115)   317     69        0 (1193/1153-1271/1224) -59.0 W12-W17 N10-N15
      7 1084 1160 (N13W05 L=106)   304     98        0 (1032/1127-1133/1186)  15.9 W01-W08 N11-N14
      8 1667 1132 (N03W44 L=144)   312    112        1 (1622/1121-1708/1148) -15.9 W40-W47 N02-N04
      9  233 1036 (N13E57 L= 43)   612     97        1 ( 202/1005- 269/1078) -33.7 E60-E53 N11-N15
     10 1536  871 (S10W31 L=131)   316     84        1 (1503/ 832-1560/ 915)  59.0 W28-W32 S12-S07
     11 1421  546 (S30W21 L=121)  1728    355        0 (1241/ 460-1565/ 667)  19.8 W07-W32 S34-S23
     12  624  475 (S25E36 L= 65)   865    220        0 ( 541/ 437- 691/ 530)  19.7 E43-E29 S27-S21
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y