Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1410 1575 (N29W43 L=157)   337     53        1 (1387/1557-1443/1594) -63.4 W40-W44 N26-N30
      2  280 1087 (N25E52 L= 62)   613    134        2 ( 255/1029- 305/1138) -74.1 E54-E49 N21-N28
      3  350 1065 (N23E45 L= 69)   360     72        0 ( 331/1029- 368/1100)  45.0 E47-E43 N20-N24
      4  382  965 (N16E43 L= 70)   437    171        1 ( 347/ 895- 439/ 997)  33.7 E46-E38 N12-N18
      5 1588  923 (S13W32 L=146)   448     96        0 (1532/ 910-1637/ 939) -15.9 W28-W35 S14-S10
      6 1767  894 (S20W47 L=161)   403     77        0 (1750/ 854-1784/ 933)  53.1 W44-W49 S21-S17
      7  920  488 (S23E20 L= 94)  3639    936        1 ( 708/ 335-1173/ 725)  18.4 E38-W01 S31-S13
      8 1299  571 (S27W06 L=120)   480    206        0 (1249/ 504-1323/ 660)  45.0 W01-W09 S30-S21
      9  684  279 (S33E48 L= 65)   650    104        1 ( 658/ 245- 709/ 321)  50.2 E53-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