Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  645 1731 (N31E48 L=236)   373     82        0 ( 618/1704- 663/1752)  -9.5 E51-E45 N29-N32
      2  982 1661 (N34E21 L=263)   344     66        0 ( 959/1642-1017/1687)   0.0 E23-E18 N32-N35
      3  415 1633 (N19E62 L=223)   490    110        1 ( 397/1602- 430/1667) -45.0 E64-E59 N17-N21
      4  746 1598 (N24E33 L=251)   414     74        0 ( 721/1568- 770/1626) -31.0 E35-E30 N22-N25
      5  361 1473 (N08E57 L=227)   944    130        3 ( 329/1418- 410/1520) -60.3 E61-E52 N04-N11
      6  265 1099 (S16E53 L=231)  1651    292        4 ( 220/1009- 309/1222)  85.9 E56-E50 S21-S07
      7  326 1135 (S12E48 L=236)   402     81        2 ( 300/1105- 359/1160)  11.3 E50-E45 S14-S10
      8  738 1116 (S05E19 L=266)   208    127        0 ( 680/1084- 797/1157)  -6.3 E23-E14 S05-S03
      9  808  923 (S14E10 L=274)   309     82        1 ( 776/ 887- 850/ 946)  36.9 E11-E07 S16-S12
     10  267  817 (S33E49 L=235)   341     64        2 ( 242/ 796- 293/ 834)   0.0 E51-E46 S34-S31
     11  766  696 (S29E08 L=276)   303     78        1 ( 731/ 660- 798/ 720)  59.0 E09-E06 S31-S26
     12  591  563 (S42E18 L=266)   642    151        2 ( 551/ 527- 629/ 601)  45.0 E21-E15 S45-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