Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1299 1797 (N58E05 L=217)   881     81        1 (1281/1749-1327/1840) -65.6 E09-E02 N52-N63
      2 1040 1620 (N33E16 L=207)   675    237        1 ( 987/1528-1109/1704) -33.7 E23-E08 N28-N38
      3 1553 1545 (N44W25 L=248)   626     91        0 (1516/1504-1590/1576) -23.2 W21-W29 N41-N46
      4 1564 1315 (N28W28 L=250)   344    122        0 (1546/1261-1593/1372) -50.2 W25-W30 N25-N31
      5 1817 1225 (N30W54 L=276)  1317    238        0 (1744/1132-1889/1303) -17.5 W44-W64 N26-N33
      6  227  769 (S37E59 L=164)   734     81        1 ( 201/ 732- 255/ 800)  32.0 E63-E53 S39-S34
      7 1461  600 (S16W39 L=262)   290    136        0 (1399/ 578-1504/ 643)  50.2 W35-W40 S19-S12
      8  287  585 (S49E57 L=166)   655     80        2 ( 264/ 553- 323/ 620)  51.3 E62-E49 S51-S46
      9 1356  437 (S29W40 L=263)  2431    480        2 (1136/ 366-1553/ 505)  35.8 W27-W52 S37-S19
  

    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