Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  258 1288 (S07E59 L= 54)   364     73        1 ( 231/1257- 290/1318)  45.0 E61-E56 S09-S05
      2 1674 1106 (N15W38 L=151)   499    130        1 (1649/1035-1693/1165) -60.3 W34-W40 N10-N17
      3 1787 1050 (N15W49 L=162)   472    169        0 (1759/1001-1845/1122) -20.0 W44-W56 N13-N18
      4  359  996 (S23E43 L= 70)   465    155        1 ( 293/ 969- 426/1024)  14.0 E49-E37 S25-S21
      5  231  776 (S39E58 L= 54)   412     70        0 ( 206/ 747- 250/ 813)  24.0 E63-E54 S40-S36
      6 1522  652 (S14W41 L=154)   527    223        0 (1444/ 573-1594/ 725)  71.6 W39-W43 S20-S08
      7 1027  611 (S31W12 L=124)   643    207        3 ( 951/ 561-1090/ 680)  68.2 W07-W14 S35-S25
  

    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