Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1336 1717 (N49E02 L=336)   778    103        2 (1299/1676-1380/1752) -29.7 E05-W02 N46-N51
      2 1068 1488 (N23E11 L=328)  1222    392        2 (1026/1324-1121/1678) -69.4 E16-E06 N12-N36
      3  704 1636 (N22E40 L=299)   556     98        1 ( 671/1600- 742/1673) -80.5 E42-E36 N19-N25
      4  542 1573 (N14E48 L=291)  1310    218        1 ( 473/1502- 605/1668) -85.6 E51-E45 N08-N21
      5 1589 1513 (N42W25 L=  4)  2171    288        2 (1471/1418-1701/1592)  -9.9 W13-W37 N39-N45
      6 1364 1466 (N31W08 L=347)   253     70        0 (1352/1426-1373/1502) -53.1 W06-W09 N28-N33
      7 1350 1326 (N21W10 L=349)   260    117        0 (1296/1293-1387/1359)  -8.1 W06-W13 N20-N21
      8 1539 1185 (N19W26 L=  5)   898    203        1 (1516/1105-1560/1289) -80.5 W23-W28 N13-N25
      9 1405 1199 (N15W17 L=356)   331    131        0 (1388/1127-1417/1255) -54.5 W14-W19 N11-N18
     10  188 1153 (S18E65 L=274)  1786    191        3 ( 152/1068- 219/1244)  84.8 E69-E61 S23-S12
     11 1210 1129 (N06W07 L=346)   248     83        0 (1181/1101-1231/1172) -76.0 W06-W08 N03-N08
     12 1750  909 (N10W49 L= 28)   433     68        3 (1729/ 872-1781/ 945) -53.1 W46-W52 N07-N12
     13  926  729 (S24W03 L=342)  1563    415        2 ( 859/ 630-1069/ 887) -32.0 E03-W14 S28-S15
     14 1162  797 (S14W14 L=353)   223     68        0 (1150/ 754-1172/ 832) -53.1 W12-W15 S15-S11
     15  543  722 (S36E23 L=316)   834    247        0 ( 439/ 697- 676/ 744)  25.3 E32-E13 S39-S30
     16 1127  662 (S22W17 L=356)   217     67        0 (1097/ 643-1168/ 678)  45.0 W15-W18 S23-S20
     17  653  576 (S42E11 L=328)  1332    313        2 ( 520/ 542- 784/ 631)  33.0 E23-E02 S47-S34
     18 1549  462 (S20W54 L= 33)   953    234        2 (1440/ 434-1652/ 485)  31.0 W46-W63 S24-S15
  

    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