Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1111 1910 (N60W49 L=229)  1087    110        3 (1073/1868-1152/1949)   9.5 W38-W68 N56-N62
      2 1437 1764 (N38W56 L=236)   450     65        0 (1403/1751-1473/1777) -14.0 W52-W60 N36-N38
      3 1912 1151 (S09W69 L=249)   537     55        1 (1899/1116-1925/1181)  26.6 W66-W72 S10-S07
      4 1790 1115 (S08W53 L=232)   352     58        1 (1772/1082-1805/1155)  76.0 W51-W54 S10-S06
      5  400 1099 (N18E39 L=141)   268     49        1 ( 375/1083- 432/1114) -11.3 E40-E35 N17-N19
      6 1293  720 (S21W10 L=190)   774    291        2 (1168/ 640-1402/ 810)  15.5 W00-W18 S23-S17
      7  838  747 (S10E16 L=163)   315    128        0 ( 822/ 680- 857/ 792)  74.1 E17-E14 S14-S07
      8  382  651 (S08E50 L=129)  1017    180        0 ( 356/ 567- 412/ 738)  73.3 E52-E47 S12-S02
      9 1129  520 (S30E04 L=175)   213     80        0 (1105/ 477-1156/ 555)  38.7 E06-E01 S32-S28
     10  973  402 (S35E18 L=161)  1022    214        3 ( 873/ 375-1090/ 436) -24.8 E25-E10 S39-S30
  

    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