Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1407 1804 (N39W62 L=288)  2268    218        5 (1315/1772-1523/1829) -13.5 W50-W75 N34-N42
      2 1660 1623 (N20W66 L=292)  2670    315        2 (1564/1524-1789/1724) -70.6 W62-W75 N11-N28
      3 1297 1622 (N31W36 L=262)   584    164        1 (1220/1603-1382/1638) -33.7 W30-W40 N27-N33
      4  973 1459 (N30W08 L=234)   496    198        2 ( 889/1427-1054/1494) -26.6 W02-W12 N25-N32
      5  629 1425 (N37E16 L=210)   949    347        0 ( 484/1374- 803/1468)  -6.6 E28-E02 N34-N38
      6 1283 1293 (N12W22 L=247)   532    135        1 (1245/1238-1330/1355) -82.9 W20-W22 N08-N16
      7  870 1262 (N21E03 L=222)   186    106        0 ( 828/1232- 910/1290)  18.4 E06-W00 N19-N21
      8 1781 1114 (S10W51 L=277)   470     92        0 (1753/1073-1806/1158) -78.7 W49-W53 S12-S07
      9  213  894 (N14E60 L=166)   970    103        2 ( 164/ 855- 259/ 933) -31.0 E65-E54 N10-N16
     10  830  861 (S01E15 L=211)   185     34        1 ( 808/ 844- 851/ 880)  18.4 E16-E13 S01-N00
     11  920  828 (S05E10 L=215)   481    147        1 ( 864/ 762- 954/ 876) -77.5 E12-E09 S09-N00
     12 1839  773 (S33W55 L=281)   778    114        2 (1819/ 730-1875/ 832) -39.3 W51-W62 S37-S28
     13 1190  638 (S22W00 L=226)   849    195        2 (1140/ 567-1226/ 731)  45.0 E04-W04 S25-S17
     14 1680  690 (S33W37 L=262)   391     69        2 (1646/ 669-1713/ 712)   0.0 W33-W40 S34-S31
     15 1580  518 (S42W27 L=253)  1481    257        1 (1465/ 444-1650/ 606)  15.5 W15-W33 S45-S37
     16  546  380 (S23E52 L=174)   482     98        0 ( 532/ 331- 566/ 432)  31.0 E56-E46 S26-S20
     17  673  328 (S30E46 L=180)   395     70        0 ( 653/ 291- 691/ 366)  26.6 E49-E41 S31-S27
     18 1153  252 (S50E17 L=209)  1681    298        1 (1025/ 200-1285/ 311)   0.0 E30-E03 S51-S47
  

    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