Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  973 1664 (N34E22 L=263)   323     45        0 ( 949/1645-1005/1681)  33.7 E24-E19 N32-N36
      2  410 1636 (N19E62 L=223)   507    110        1 ( 392/1602- 426/1667) -51.3 E65-E60 N16-N21
      3  745 1560 (N22E32 L=253)   771    162        1 ( 713/1491- 774/1630) -41.6 E36-E27 N17-N25
      4  471 1594 (N18E54 L=231)   387     72        0 ( 453/1561- 494/1629) -78.7 E54-E52 N15-N20
      5  352 1480 (N08E59 L=226)  1032    139        3 ( 320/1439- 406/1527) -56.3 E62-E53 N05-N11
      6  258 1124 (S14E54 L=231)  1351    166        3 ( 222/1051- 294/1213) -84.3 E58-E51 S18-S08
      7  315 1141 (S12E49 L=236)   390     74        2 ( 289/1111- 345/1167) -76.0 E51-E46 S14-S10
      8  192 1068 (S19E60 L=225)   876    131        2 ( 165/ 998- 215/1112)  81.9 E62-E57 S23-S16
      9  798  934 (S14E11 L=274)   354     71        1 ( 765/ 901- 834/ 971)  78.7 E12-E08 S16-S11
     10  262  819 (S33E50 L=235)   418     72        2 ( 236/ 795- 289/ 838)   0.0 E52-E47 S34-S31
     11  756  704 (S28E09 L=277)   337     96        0 ( 716/ 668- 793/ 741)  63.4 E10-E06 S31-S25
     12  582  570 (S42E19 L=266)   594    100        1 ( 537/ 531- 617/ 607)  40.6 E23-E16 S45-S38
  

    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