Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1116 1913 (N60W50 L=229)  1067     97        2 (1083/1872-1153/1949)   9.8 W39-W68 N56-N62
      2 1440 1764 (N38W56 L=235)   432     48        0 (1410/1753-1471/1777)  -8.1 W52-W59 N36-N38
      3  381 1181 (N23E40 L=139)   253     50        0 ( 359/1159- 401/1207) -53.1 E41-E37 N21-N25
      4 1915 1148 (S09W69 L=248)   537     54        1 (1898/1116-1929/1178)  26.6 W66-W73 S10-S07
      5 1795 1115 (S08W53 L=232)   368     48        0 (1776/1084-1811/1149) -90.0 W51-W54 S10-S06
      6  406 1102 (N18E38 L=141)   270     43        0 ( 381/1087- 438/1116)   0.0 E40-E35 N17-N19
      7 1295  719 (S21W10 L=189)   708    292        2 (1175/ 645-1405/ 810)  13.2 W01-W18 S23-S17
      8  845  750 (S10E16 L=163)   320    107        0 ( 829/ 687- 864/ 794)  74.1 E17-E13 S14-S07
      9  386  654 (S07E50 L=129)  1018    179        0 ( 359/ 567- 417/ 740)  68.2 E51-E47 S12-S02
     10 1139  531 (S30E03 L=175)   233     84        0 (1112/ 489-1165/ 563)  45.0 E06-E01 S32-S27
     11  978  403 (S35E18 L=161)   986    217        2 ( 876/ 376-1092/ 438) -26.6 E24-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