Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1626 1540 (N37W36 L=132)   780    117        2 (1572/1512-1668/1560)  12.5 W31-W40 N35-N39
      2  365 1337 (N05E48 L= 48)   515    134        0 ( 305/1319- 412/1357)  12.5 E53-E44 N03-N06
      3  200 1179 (S06E61 L= 35)   446    112        2 ( 180/1138- 230/1218) -45.0 E64-E57 S08-S03
      4 1331 1158 (N06W16 L=112)   423    139        2 (1251/1119-1389/1190)  32.0 W11-W19 N02-N08
      5  224 1081 (S12E57 L= 39)  1032    257        4 ( 174/1014- 273/1184)  51.3 E61-E52 S15-S05
      6  136 1103 (S11E68 L= 28)   822    286        0 (  84/1083- 175/1126)   0.0 E77-E62 S11-S09
      7 1652 1003 (N03W40 L=136)   889    197        2 (1619/ 899-1677/1086)  84.8 W37-W41 S03-N08
      8  301  963 (S18E48 L= 48)  1806    532        5 ( 173/ 865- 421/1082)  33.7 E60-E39 S24-S10
      9 1256  918 (S09W15 L=111)   914    413        3 (1148/ 762-1363/1079)  70.7 W11-W19 S19-N01
     10 1855 1015 (N08W58 L=154)   558    107        1 (1823/ 981-1878/1047)  59.0 W55-W60 N05-N10
     11 1086  928 (S11W05 L=101)   644    269        0 ( 956/ 875-1197/ 968)  33.7 E01-W11 S15-S07
     12  542  870 (S22E29 L= 67)   424    142        2 ( 487/ 830- 597/ 923)  40.6 E32-E25 S24-S17
     13 1308  826 (S13W20 L=116)   249     95        2 (1284/ 774-1343/ 850) -59.0 W18-W22 S16-S11
     14  400  723 (S32E40 L= 56)   286     63        0 ( 370/ 706- 429/ 739)   0.0 E42-E37 S33-S31
     15 1569  468 (S29W50 L=146)  1754    255        2 (1448/ 405-1684/ 548)  47.3 W43-W56 S35-S21
     16 1222  311 (S47W33 L=129)   319     67        0 (1195/ 289-1242/ 342)  78.7 W31-W34 S49-S44
  

    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