Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  452 1665 (N37E61 L=302)  2021    148        4 ( 422/1598- 487/1737) -16.2 E82-E50 N32-N41
      2  651 1614 (N36E36 L=327)  1948    231        4 ( 551/1576- 748/1657)   0.0 E44-E26 N31-N38
      3  615 1405 (N21E31 L=331)   664    219        1 ( 548/1357- 659/1515)  74.7 E35-E29 N17-N28
      4 1222 1387 (N26W10 L= 12)   341    142        0 (1195/1333-1277/1421)  56.3 W08-W13 N22-N28
      5 1714 1302 (N24W49 L= 51)   306     73        1 (1698/1261-1729/1336) -78.7 W47-W50 N21-N26
      6 1014 1320 (N20E03 L=359)   274    116        1 ( 959/1301-1076/1330)   9.5 E06-W00 N18-N20
      7 1439 1194 (N15W25 L= 28)   208     78        0 (1399/1184-1479/1213)  -9.5 W22-W28 N14-N16
      8 1770 1070 (N10W52 L= 55)   589    148        0 (1746/ 987-1792/1120) -63.4 W49-W54 N05-N13
      9 1430  995 (N03W25 L= 28)   473    144        1 (1376/ 951-1496/1027)  20.6 W22-W30 N00-N05
     10 1758  721 (S11W56 L= 59)  3471    503        3 (1626/ 586-1870/ 943)  54.0 W47-W64 S20-N02
     11  169  866 (S16E66 L=296)   653     82        1 ( 137/ 837- 198/ 903) -56.3 E72-E61 S18-S14
     12 1529  874 (S04W34 L= 36)   256     49        1 (1510/ 852-1557/ 897)  18.4 W32-W35 S04-S02
     13  433  594 (S32E41 L=322)  1109    226        0 ( 364/ 512- 517/ 648)  31.0 E52-E32 S38-S26
  

    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