Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  701 1696 (N38E31 L=354)  2055    307        4 ( 553/1637- 828/1736) -12.1 E46-E18 N33-N40
      2 1222 1676 (N40W10 L= 35)   731    214        0 (1123/1637-1312/1730) -21.8 W02-W17 N36-N43
      3 1052 1379 (N17E01 L= 25)  1029    246        1 ( 990/1277-1098/1485)  77.9 E03-W02 N10-N24
      4  379 1423 (N17E48 L=337)   308     91        0 ( 346/1412- 432/1443) -12.5 E52-E43 N16-N18
      5  257 1335 (N12E58 L=327)   384     54        0 ( 238/1309- 279/1368) -33.7 E61-E54 N09-N13
      6 1256 1290 (N13W13 L= 38)   551    242        0 (1149/1241-1362/1338) -12.1 W05-W19 N10-N15
      7  471 1279 (N08E37 L=348)   296     92        1 ( 432/1255- 518/1294)  15.9 E40-E33 N06-N09
      8  276 1168 (N01E52 L=333)   396     78        0 ( 248/1132- 305/1193)  21.8 E55-E49 N00-N02
      9  522 1111 (S02E32 L=353)   423    116        0 ( 480/1066- 568/1156)  45.0 E34-E29 S04-N01
     10 1234 1081 (N00W12 L= 37)   209     71        0 (1198/1062-1265/1099)   0.0 W10-W14 N00-N01
     11  254 1007 (S08E54 L=331)   653    116        2 ( 237/ 943- 281/1048)  60.3 E55-E51 S12-S05
     12  179  961 (S11E63 L=322)   569     79        3 ( 161/ 922- 198/1002)  78.7 E65-E60 S13-S08
     13  136  922 (S13E70 L=315)   660     99        2 ( 121/ 883- 163/ 981)  63.4 E73-E65 S15-S09
     14 1536  864 (S10W34 L= 59)   278    116        0 (1525/ 793-1550/ 905)  90.0 W33-W34 S14-S07
     15  946  851 (S15E04 L= 21)   390    202        1 ( 842/ 829-1029/ 871)   5.2 E10-W01 S16-S14
     16 1777  671 (S19W60 L= 86)  3235    412        9 (1700/ 561-1852/ 778)  65.0 W53-W70 S26-S11
     17  216  637 (S31E69 L=316)  2507    185        2 ( 188/ 541- 240/ 722)  32.9 E83-E63 S36-S25
     18 1632  620 (S24W48 L= 73)   608    178        0 (1585/ 543-1663/ 687) -36.0 W41-W52 S28-S20
  

    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