Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1062 1627 (N34E01 L=259)   278     44        0 (1040/1601-1092/1645)  56.3 E02-W01 N31-N35
      2 1657 1577 (N35W48 L=308)   485     65        0 (1626/1554-1690/1595)  12.5 W43-W52 N33-N36
      3 1554 1543 (N32W36 L=296)   416     75        0 (1541/1503-1569/1585)  56.3 W34-W38 N28-N34
      4  284 1483 (N22E59 L=200)  1093    101        4 ( 258/1436- 321/1532) -66.8 E63-E55 N18-N25
      5 1834 1454 (N29W65 L=325)   837    105        1 (1800/1426-1872/1479)  -5.7 W59-W72 N26-N30
      6  261 1370 (N14E57 L=203)   486     88        0 ( 204/1338- 294/1396)  16.7 E63-E53 N12-N16
      7  441 1341 (N13E40 L=220)   347    136        1 ( 396/1292- 505/1366)  35.5 E42-E35 N09-N14
      8  640 1274 (N09E25 L=235)   435    147        0 ( 573/1256- 713/1298)   0.0 E29-E20 N07-N10
      9   94  974 (S09E75 L=185)  1035     92        1 (  77/ 929- 120/1009)  39.8 E79-E69 S11-S06
     10  173  833 (S18E64 L=196)   457     71        2 ( 160/ 786- 186/ 869)  59.0 E65-E61 S20-S15
     11 1527  597 (S27W38 L=298)  1633    279        4 (1434/ 494-1641/ 673)  53.1 W32-W45 S34-S21
     12 1335  579 (S30W24 L=284)   386    122        0 (1284/ 550-1383/ 622)  45.0 W21-W26 S31-S26
     13  671  569 (S34E23 L=237)   862    272        0 ( 554/ 483- 752/ 613) -23.2 E32-E17 S39-S31
     14 1651  522 (S30W53 L=312)   472     84        0 (1609/ 492-1686/ 559)  59.0 W49-W55 S32-S27
     15  452  516 (S38E44 L=216)   372     78        2 ( 420/ 489- 495/ 545)  16.7 E48-E38 S39-S35
     16 1606  327 (S44W66 L=325)  1027    108        1 (1562/ 298-1654/ 364)  -2.2 W56-W83 S44-S41
     17 1016  317 (S51W05 L=264)  1411    286        2 ( 863/ 294-1131/ 352)   4.8 E09-W15 S53-S48
  

    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