Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  976 1663 (N35E22 L=263)   349     51        0 ( 951/1645-1010/1687)   0.0 E24-E19 N32-N35
      2  749 1568 (N22E32 L=253)   666    164        2 ( 716/1497- 780/1629) -37.9 E36-E27 N18-N25
      3  473 1589 (N17E53 L=232)   357     64        0 ( 454/1559- 497/1618)  90.0 E54-E51 N15-N19
      4  357 1480 (N08E58 L=226)   935    177        4 ( 327/1433- 409/1526) -45.0 E62-E52 N05-N11
      5  605 1424 (N10E36 L=249)   242     51        0 ( 582/1401- 629/1451)  76.0 E36-E34 N08-N12
      6  261 1101 (S16E53 L=231)  1493    269        4 ( 223/1014- 301/1216)  77.0 E57-E50 S21-S08
      7  324 1134 (S12E48 L=236)   407     96        1 ( 298/1102- 357/1163) -53.1 E50-E45 S14-S10
      8  197 1069 (S19E60 L=225)   815    131        3 ( 179/1011- 218/1113) -80.5 E62-E57 S22-S16
      9  802  927 (S14E11 L=274)   310     70        0 ( 771/ 898- 841/ 947)  45.0 E12-E08 S16-S12
     10  265  817 (S33E49 L=235)   368     68        2 ( 240/ 797- 291/ 834)   9.5 E52-E46 S34-S31
     11  764  703 (S28E08 L=277)   334    102        0 ( 722/ 667- 796/ 740)  56.3 E10-E06 S31-S25
     12  587  566 (S42E19 L=266)   654    129        1 ( 538/ 532- 628/ 604)  36.9 E23-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