Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  542 1734 (N70E28 L=110)   597     59        3 ( 519/1700- 566/1764)  74.1 E35-E23 N66-N73
      2  475 1704 (N70E43 L= 95)   508     47        0 ( 452/1686- 495/1722)   0.0 E50-E36 N67-N71
      3 1110 1651 (N41W28 L=166)   259     43        2 (1078/1633-1135/1669) -11.3 W25-W30 N39-N42
      4 1499 1555 (N23W48 L=186)   286     80        0 (1465/1523-1531/1580) -68.2 W47-W49 N20-N25
      5 1666 1512 (N15W60 L=198)   615     71        2 (1643/1475-1691/1539)  14.0 W56-W64 N13-N16
      6 1578 1496 (N17W51 L=189)   286     75        1 (1544/1468-1603/1516) -21.8 W48-W53 N15-N18
      7 1290 1383 (N20W26 L=164)   465    311        0 (1148/1329-1406/1418) -40.2 W18-W31 N13-N25
      8  730 1287 (N29E11 L=127)   271    118        0 ( 694/1247- 782/1328) -56.3 E13-E08 N26-N32
      9  246 1287 (N40E57 L= 81)   466     40        1 ( 223/1263- 275/1308) -14.0 E61-E53 N38-N42
     10  608 1194 (N27E22 L=116)   336    111        1 ( 575/1161- 656/1237) -50.2 E23-E18 N24-N30
     11  214  950 (N20E61 L= 77)   366     64        0 ( 187/ 922- 241/ 971)  15.9 E64-E57 N18-N21
     12  383  881 (N13E45 L= 93)   643    156        1 ( 323/ 858- 455/ 914) -31.0 E49-E38 N10-N16
     13  244  882 (N15E58 L= 80)   402     77        0 ( 220/ 850- 262/ 913)  45.0 E61-E55 N13-N17
     14  596  835 (N06E30 L=108)   445    122        2 ( 569/ 773- 643/ 880)  74.1 E32-E26 N02-N09
     15  890  819 (S01E13 L=125)   253     70        0 ( 866/ 790- 913/ 848) -71.6 E13-E11 S03-N00
     16  558  540 (S10E42 L= 96)   415     65        0 ( 532/ 505- 577/ 578)  78.7 E43-E40 S12-S07
     17 1375  509 (S34W08 L=146)   642     73        2 (1334/ 477-1411/ 543)  -9.5 W05-W11 S36-S30
     18  387  454 (S13E63 L= 75)  2050    175        7 ( 323/ 395- 454/ 514) -52.1 E72-E56 S17-S08
  

    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