Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1160 1644 (N41W28 L=338)   266     71        0 (1141/1615-1179/1670)  53.1 W26-W30 N38-N42
      2 1219 1500 (N30W25 L=335)   272     64        1 (1196/1478-1259/1521)   0.0 W22-W27 N29-N31
      3 1771 1442 (N11W65 L= 15)  1592    198        2 (1721/1385-1821/1510) -77.5 W61-W71 N07-N16
      4 1532 1211 (N06W35 L=345)   744    290        0 (1510/1070-1567/1339)  71.6 W31-W38 S02-N13
      5  148  992 (N19E66 L=244)  1491    320        2 (  93/ 936- 202/1051) -32.5 E72-E60 N15-N22
      6 1653  814 (S20W36 L=346)   260     22        0 (1637/ 793-1668/ 836)  33.7 W34-W37 S21-S18
      7 1392  458 (S37W12 L=322)   306     67        0 (1371/ 439-1420/ 484) -59.0 W10-W13 S39-S34
  

    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