Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  834 1676 (N34E30 L=181)  2676    565        9 ( 635/1569-1027/1777)  -6.2 E50-E13 N30-N38
      2 1109 1555 (N32E07 L=204)   464     90        0 (1074/1521-1166/1577)  50.2 E09-E03 N29-N35
      3  580 1428 (N13E37 L=174)   313     71        0 ( 547/1401- 603/1465)  68.2 E38-E35 N10-N15
      4  399 1398 (N08E50 L=161)   976    276        2 ( 302/1362- 504/1441)  16.4 E59-E42 N04-N10
      5  242 1324 (N01E63 L=147)   622    120        0 ( 218/1300- 275/1349)  -6.3 E67-E58 N00-N01
      6 1758 1026 (N14W49 L=260)  1033    210        1 (1691/ 951-1801/1085) -24.4 W42-W54 N10-N16
      7  725 1008 (S09E17 L=193)   533    150        0 ( 642/ 981- 792/1042)  29.1 E22-E13 S11-S05
      8  244 1010 (S18E54 L=157)   417    111        1 ( 213/ 991- 278/1027)  18.4 E57-E50 S19-S16
      9 1513  714 (S10W38 L=248)   296     90        0 (1497/ 675-1537/ 759) -45.0 W35-W40 S11-S07
     10  687  671 (S30E14 L=197)   276     95        1 ( 634/ 652- 732/ 690)   7.1 E18-E10 S31-S28
     11  948  639 (S26W04 L=215)   446     90        1 ( 896/ 619- 996/ 670)  15.9 W00-W07 S27-S24
  

    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