Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  574 1621 (N53E31 L=240)   361     51        0 ( 561/1592- 592/1663) -60.3 E34-E28 N50-N57
      2  657 1621 (N52E21 L=250)   412     47        1 ( 632/1590- 686/1642)  36.9 E24-E17 N49-N53
      3 1793 1369 (N12W65 L=336)   840    134        2 (1772/1316-1818/1430)  81.9 W61-W69 N08-N15
      4 1170 1314 (N21W14 L=285)   337    156        1 (1104/1283-1241/1355) -40.6 W10-W18 N18-N24
      5  965  848 (S03E06 L=265)   199     93        0 ( 916/ 840-1010/ 859)  -9.5 E09-E03 S04-S02
      6 1147  475 (S30E01 L=270)  2680    553        2 ( 938/ 346-1306/ 656)  28.1 E19-W11 S37-S20
      7  437  363 (S29E64 L=207)   986    134        0 ( 417/ 305- 457/ 430)  26.6 E78-E55 S35-S24
      8  595  376 (S29E45 L=226)   537    106        0 ( 577/ 329- 616/ 428)  63.4 E47-E42 S33-S25
  

    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