Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1111 1715 (N39E18 L= 97)  2386    476        5 ( 960/1580-1207/1825) -17.7 E33-E08 N29-N43
      2  551 1644 (N18E50 L= 64)  2366    438       11 ( 467/1561- 623/1737)  90.0 E56-E43 N12-N24
      3  347 1531 (N07E62 L= 52)  1400    197        2 ( 303/1467- 391/1609) -60.9 E67-E58 N02-N11
      4  686 1484 (N11E33 L= 82)   879    270        1 ( 612/1395- 732/1544) -22.6 E39-E27 N07-N15
      5  763 1458 (N12E27 L= 87)   392    141        0 ( 740/1391- 804/1515) -77.5 E28-E25 N07-N16
      6 1545 1225 (N21W25 L=140)   778    125        1 (1491/1182-1587/1283) -38.7 W22-W28 N16-N23
      7 1220  237 (S45W47 L=162)  1620    268        2 (1101/ 212-1335/ 273)  26.6 W38-W55 S48-S39
  

    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