Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1398 1570 (N29W42 L=157)   345     49        1 (1377/1544-1419/1596)  18.4 W38-W44 N27-N30
      2  281 1069 (N24E52 L= 63)   409     53        1 ( 263/1030- 302/1101)  63.4 E54-E49 N21-N25
      3  339 1040 (N21E46 L= 69)   468    142        1 ( 320/ 974- 375/1089) -81.9 E48-E43 N17-N24
      4 1568  916 (S12W31 L=145)   416    114        1 (1506/ 900-1619/ 930) -14.0 W26-W34 S14-S10
      5 1756  890 (S20W46 L=161)   398     57        0 (1738/ 853-1774/ 922)  45.0 W43-W47 S21-S17
      6 1038  554 (S21E10 L=104)  2042    473        3 ( 869/ 412-1158/ 720)  27.5 E25-W00 S27-S13
      7 1281  559 (S27W05 L=120)   639    221        1 (1225/ 486-1309/ 657)  51.3 W00-W08 S30-S20
      8  756  404 (S25E34 L= 81)  1701    436        0 ( 693/ 331- 875/ 551) -82.4 E39-E26 S31-S15
      9  675  276 (S33E49 L= 66)   659    105        1 ( 650/ 243- 701/ 319)  50.2 E54-E45 S35-S29
  

    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