Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1892 (N52E35 L=292)   628     74        1 ( 890/1876- 964/1914) -10.3 E42-E30 N49-N54
      2  680 1792 (N37E46 L=281)   834     95        3 ( 655/1747- 712/1846) -48.4 E53-E41 N33-N42
      3 1271 1435 (N22W08 L=335)   300     94        1 (1243/1398-1308/1459) -31.0 W05-W10 N20-N23
      4  457  929 (S21E35 L=292)   598    178        1 ( 400/ 870- 505/ 986)  48.8 E39-E31 S24-S16
      5 1713  641 (S14W56 L= 23)   483     43        2 (1688/ 617-1739/ 670)  90.0 W53-W59 S15-S12
      6  357  578 (S45E49 L=279)   712     77        1 ( 336/ 540- 377/ 624)  80.5 E52-E45 S47-S41
      7 1694  521 (S21W62 L= 29)  2317    308        8 (1608/ 437-1779/ 599)  65.2 W55-W67 S27-S14
      8 1151  526 (S34W19 L=347)   449    153        0 (1072/ 510-1228/ 546)  33.7 W14-W23 S36-S30
      9 1025  244 (S57W26 L=353)   919    153        1 ( 947/ 215-1088/ 280)   3.4 W15-W32 S59-S54
  

    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