Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1212 1780 (N48W21 L=281)  1318    209        6 (1143/1731-1297/1841) -36.0 W15-W27 N43-N54
      2 1610 1660 (N38W54 L=315)  1246    183        2 (1515/1639-1691/1681)   0.0 W43-W70 N36-N39
      3 1542 1493 (N25W38 L=299)  6329    971        9 (1144/1240-1878/1623) -21.4 W10-W66 N09-N35
      4 1524 1369 (N17W34 L=295)   458    112        0 (1508/1313-1549/1411) -90.0 W32-W36 N13-N20
      5  816 1345 (N17E12 L=249)   561    123        3 ( 765/1302- 866/1382) -45.0 E15-E08 N14-N20
      6 1743 1249 (N10W50 L=311)   432     82        1 (1728/1200-1769/1291)  71.6 W48-W53 N06-N12
      7 1653 1258 (N10W42 L=303)   287     49        0 (1629/1236-1677/1287)  53.1 W40-W44 N08-N12
      8 1833 1217 (N08W59 L=320)   718     60        1 (1807/1186-1860/1247) -36.9 W56-W62 N06-N09
      9  403 1204 (N11E40 L=220)   348    130        0 ( 354/1183- 469/1234)  12.5 E44-E35 N09-N12
     10  967 1193 (N08E03 L=258)   276     92        0 ( 928/1157-1003/1229) -45.0 E05-W00 N05-N09
     11 1886  875 (S13W66 L=326)  2832    283        1 (1826/ 759-1945/1012)  53.7 W59-W74 S19-S04
     12 1600  797 (S18W38 L=299)   844    276        0 (1558/ 685-1663/ 925)  71.6 W35-W42 S24-S09
     13  404  766 (S16E43 L=218)   831    211        4 ( 306/ 739- 498/ 791)  -3.6 E51-E35 S17-S13
     14  389  601 (S26E49 L=211)   581    114        0 ( 342/ 568- 436/ 643)  18.4 E55-E43 S27-S23
     15 1266  590 (S30W15 L=276)   226     74        0 (1240/ 557-1301/ 621)  53.1 W13-W17 S32-S28
     16 1734  506 (S36W62 L=323)  1716    195        1 (1649/ 472-1806/ 560)  16.4 W50-W72 S38-S32
     17 1467  464 (S39W34 L=295)   683    109        3 (1417/ 422-1510/ 505) -12.5 W29-W38 S42-S36
     18 1149  449 (S40W08 L=268)   769    179        2 (1066/ 419-1233/ 486)  14.0 W01-W13 S41-S37
     19  374  384 (S40E69 L=192)   647     72        0 ( 342/ 359- 397/ 416)   7.8 E82-E60 S41-S38
  

    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