Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  675 1699 (N59E19 L= 14)  1263    191        3 ( 611/1646- 764/1744)   9.0 E27-E08 N54-N62
      2  550 1648 (N57E35 L=359)  1041    133        2 ( 511/1597- 602/1690)  82.9 E41-E28 N52-N60
      3  995 1604 (N44W11 L= 44)   717    148        1 ( 922/1580-1073/1629) -31.0 W04-W16 N40-N47
      4 1216 1585 (N38W27 L= 60)   676    129        2 (1157/1561-1287/1603) -11.3 W22-W32 N36-N39
      5 1216 1370 (N24W19 L= 52)   253    149        0 (1163/1353-1287/1400) -29.7 W15-W23 N22-N26
      6  335 1291 (N33E49 L=345)   320     48        1 ( 310/1271- 364/1311)   9.5 E51-E45 N31-N33
      7  956 1187 (N17E01 L= 32)   369    136        0 ( 908/1136- 988/1233) -66.8 E03-W00 N13-N20
      8  442 1181 (N24E38 L=356)   308     75        1 ( 427/1147- 460/1226) -78.7 E38-E36 N22-N27
      9  316 1035 (N17E49 L=344)   284     32        0 ( 294/1019- 341/1051)  11.3 E51-E46 N15-N17
     10 1788  730 (S27W53 L= 86)  2047    397        2 (1671/ 622-1912/ 895)  20.8 W40-W70 S32-S20
     11 1359  650 (S22W15 L= 48)  1630    291        6 (1293/ 540-1432/ 784) -77.5 W10-W19 S31-S13
     12  723  546 (S18E28 L=  5)   287     76        0 ( 700/ 511- 741/ 581)  45.0 E30-E25 S20-S16
  

    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