Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  597 1739 (N30E50 L= 65)  1338    201        2 ( 521/1673- 671/1794) -15.6 E64-E39 N26-N34
      2  949 1686 (N33E20 L= 95)   796    128        2 ( 911/1624- 981/1745) -52.1 E24-E16 N28-N37
      3  851 1590 (N24E23 L= 92)   230     94        0 ( 809/1562- 887/1611)  53.1 E24-E21 N21-N25
      4 1072 1484 (N21E07 L=109)   291     59        2 (1049/1455-1107/1510) -14.0 E08-E04 N19-N22
      5 1708 1395 (N31W41 L=156)   364    102        0 (1669/1348-1730/1425) -38.7 W36-W42 N28-N33
      6 1614 1366 (N27W32 L=147)   533     85        1 (1592/1323-1640/1414) -59.0 W29-W34 N24-N29
      7 1211 1374 (N17W04 L=119)   273     63        1 (1182/1348-1248/1401)  53.1 W02-W06 N15-N19
      8  884 1281 (N05E13 L=102)   737    246        2 ( 861/1173- 911/1393) -77.0 E14-E09 S01-N12
      9 1263 1346 (N17W08 L=123)   205     41        0 (1234/1326-1283/1362)   0.0 W05-W09 N15-N17
     10 1846 1151 (N20W55 L=171)   539     99        1 (1802/1115-1878/1184) -12.5 W50-W59 N18-N20
     11 1409  926 (S05W25 L=140)   269     60        0 (1384/ 902-1435/ 948)  56.3 W23-W26 S06-S03
     12 1223  855 (S13W15 L=130)   340    160        1 (1152/ 824-1296/ 888)  40.6 W11-W18 S15-S09
     13 1042  805 (S19W05 L=121)   251     55        1 (1023/ 773-1068/ 831) -45.0 W03-W07 S20-S17
     14  756  792 (S25E12 L=103)   313    103        0 ( 708/ 760- 801/ 822) -15.9 E15-E08 S26-S24
     15  695  653 (S35E14 L=101)   289     74        1 ( 668/ 625- 733/ 676)  31.0 E16-E11 S37-S33
     16 1712  649 (S12W56 L=171)   560     85        0 (1679/ 627-1747/ 672)  26.6 W52-W59 S14-S10
     17 1458  469 (S29W45 L=160)  1126    297        2 (1324/ 405-1557/ 505)   2.4 W33-W57 S32-S26
     18  807  462 (S46E02 L=113)   365     49        2 ( 780/ 436- 829/ 491) -53.1 E04-W00 S47-S43
     19  830  214 (S67W11 L=126)   372     63        1 ( 794/ 193- 851/ 232)  -5.2 W05-W16 S68-S65
  

    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