Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  758 1620 (N38E27 L=309)  5603    791       10 ( 464/1519-1081/1770)  -5.7 E80-W00 N32-N45
      2 1205 1431 (N29W09 L=345)   405    121        0 (1171/1386-1249/1476) -39.8 W06-W12 N25-N31
      3  261 1428 (N20E66 L=270)   667     48        0 ( 239/1401- 291/1451)  36.9 E68-E62 N17-N21
      4  904 1364 (N22E11 L=326)   550    160        0 ( 877/1299- 943/1450) -59.0 E13-E07 N17-N27
      5 1620 1376 (N28W41 L= 17)   315     55        0 (1590/1354-1645/1393) -21.8 W38-W43 N26-N29
      6  985 1328 (N20E05 L=331)   222     99        1 ( 952/1294-1012/1370)  68.2 E06-E03 N17-N22
      7  180 1211 (N06E67 L=269)  3462    234        5 ( 147/1099- 217/1314) -67.4 E73-E60 S01-N11
      8 1806 1234 (N20W58 L= 34)   502     87        1 (1775/1192-1821/1271)  78.7 W54-W59 N17-N22
      9 1770 1161 (N15W53 L= 29)   451     83        0 (1730/1142-1811/1176)   0.0 W48-W57 N13-N16
     10  995 1057 (N03E02 L=334)   191     98        0 ( 956/1029-1035/1089) -36.9 E04-W00 N02-N05
     11 1778  961 (N03W54 L= 30)   565    129        2 (1736/ 905-1827/ 983)  24.0 W49-W58 N00-N04
     12 1380  937 (S01W23 L=359)   234     85        0 (1349/ 919-1422/ 959) -11.3 W20-W25 S02-N00
     13  427  863 (S13E39 L=298)   289     94        1 ( 392/ 821- 464/ 890)  26.6 E41-E35 S15-S11
     14  354  821 (S16E45 L=291)   342    100        0 ( 317/ 795- 393/ 858) -26.6 E48-E41 S17-S14
     15 1843  828 (S05W63 L= 39)   448     38        0 (1811/ 810-1867/ 842)   8.1 W59-W66 S05-S03
     16 1127  706 (S17W09 L=345)  1128    401        1 ( 979/ 648-1320/ 801)  31.0 W00-W20 S21-S09
     17  782  546 (S31E13 L=323)   991    190        5 ( 690/ 504- 866/ 576)  28.3 E20-E06 S34-S27
     18  520  404 (S45E39 L=297)   504    100        1 ( 486/ 384- 570/ 432) -15.9 E42-E34 S46-S42
  

    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