Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  604 1617 (N28E44 L= 12)  1001    331        3 ( 575/1466- 667/1738) -39.0 E54-E32 N18-N36
      2 1526 1495 (N39W28 L= 85)   338     70        2 (1496/1462-1556/1519)  45.0 W25-W31 N36-N41
      3 1319 1437 (N30W12 L= 68)   446    218        0 (1235/1400-1408/1505) -12.1 W04-W18 N28-N33
      4  771 1441 (N20E24 L= 32)   491    112        3 ( 724/1415- 826/1468)  26.6 E27-E21 N17-N21
      5 1671 1252 (N25W41 L= 97)   280     48        0 (1651/1228-1688/1271) -56.3 W39-W42 N23-N26
      6  198 1224 (S03E66 L=351)   687     98        3 ( 175/1174- 224/1267) -32.0 E71-E61 S06-S01
      7 1706 1155 (N20W44 L=101)   555     78        1 (1663/1133-1747/1180)   0.0 W40-W48 N18-N21
      8  425  857 (S22E37 L= 19)  2990    690        5 ( 207/ 707- 664/1091)  37.3 E65-E22 S35-S03
      9 1759 1064 (N15W50 L=107)   452     69        1 (1724/1042-1788/1082)  23.2 W46-W53 N13-N16
     10  176  968 (S19E65 L=352)  2597    294        6 ( 124/ 881- 210/1075)  40.9 E76-E60 S25-S12
     11  132  987 (S19E73 L=343)  1021    110        1 ( 115/ 944- 151/1043) -50.2 E77-E69 S21-S15
     12 1710  931 (N06W47 L=104)   321    100        0 (1682/ 887-1741/ 968) -23.2 W44-W51 N04-N07
     13 1519  883 (S00W33 L= 90)   519    110        1 (1466/ 855-1562/ 908)  18.4 W30-W36 S02-N01
     14 1514  741 (S09W37 L= 93)   279    169        0 (1483/ 715-1568/ 798) -38.7 W34-W41 S10-S05
     15 1249  705 (S16W20 L= 77)   360    137        1 (1189/ 684-1315/ 720)  35.5 W16-W23 S18-S13
     16 1716  585 (S14W61 L=118)  1628    288        1 (1625/ 500-1803/ 682)  57.3 W55-W66 S20-S06
     17  681  531 (S39E15 L= 42)   300     63        2 ( 665/ 503- 698/ 557) -56.3 E16-E13 S40-S36
     18 1492  395 (S31W52 L=108)   852     80        2 (1449/ 365-1539/ 430)  33.7 W47-W55 S33-S28
  

    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