Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  739 1734 (N34E42 L= 94)   399     74        0 ( 716/1708- 762/1764) -78.7 E43-E39 N31-N36
      2 1178 1682 (N42E09 L=126)  1119    192        1 (1134/1600-1233/1745) -31.6 E15-E02 N37-N47
      3  774 1498 (N19E28 L=108)   696    188        3 ( 730/1406- 831/1594) -81.3 E29-E25 N12-N25
      4  488 1493 (N12E47 L= 88)  1071    331        5 ( 450/1395- 538/1579) -67.4 E50-E45 N05-N17
      5  161 1251 (S08E72 L= 64)  1685    259        7 ( 126/1202- 195/1290)  45.0 E78-E66 S11-S05
      6  660 1169 (S03E25 L=111)   271    113        1 ( 631/1133- 680/1217) -78.7 E27-E22 S05-N00
      7  934  864 (S14E01 L=134)   199     50        2 ( 905/ 843- 957/ 886)  56.3 E02-W00 S15-S12
      8 1536  512 (S19W48 L=184)   387    154        4 (1508/ 474-1568/ 554) -45.0 W45-W50 S21-S17
      9 1282  326 (S38W41 L=177)  1842    369        4 (1093/ 292-1482/ 370)  26.6 W27-W53 S42-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