Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  725 1608 (N29E27 L=126)   539    163        1 ( 694/1552- 762/1672) -53.1 E28-E22 N24-N33
      2  807 1577 (N27E20 L=133)   235     69        0 ( 778/1555- 851/1594)  36.9 E21-E17 N25-N28
      3 1744 1386 (N25W49 L=201)  1171    171        2 (1670/1328-1798/1449) -26.6 W42-W53 N22-N28
      4  335 1366 (N11E50 L=103)   479    103        0 ( 308/1320- 370/1416)  63.4 E52-E47 N07-N14
      5 1254 1346 (N16W11 L=164)   254     54        0 (1218/1326-1281/1362) -14.0 W08-W12 N14-N17
      6 1358 1323 (N16W18 L=171)   259     72        0 (1326/1301-1382/1357)  63.4 W15-W19 N14-N18
      7 1142 1221 (N07W05 L=158)  1243    364        1 (1000/1086-1223/1342) -42.9 E04-W10 N00-N13
      8 1685 1312 (N19W42 L=195)   280     45        0 (1669/1287-1704/1341) -76.0 W40-W43 N17-N21
      9  581 1134 (S02E28 L=125)   624    158        1 ( 524/1085- 652/1171) -18.4 E31-E22 S04-N00
     10 1046 1046 (S04W01 L=154)   193     63        1 (1019/1030-1080/1065) -33.7 W00-W03 S04-S02
     11  891  968 (S10E07 L=145)   227     80        0 ( 852/ 947- 926/ 998)  45.0 E09-E05 S11-S07
     12  148  885 (S18E67 L= 86)  1080    133        1 ( 119/ 825- 166/ 938)  66.8 E72-E63 S21-S14
     13 1039  771 (S20W03 L=156)   998    416        0 ( 878/ 709-1207/ 850)  21.3 E06-W12 S24-S13
     14 1380  759 (S17W25 L=178)   283    127        0 (1312/ 748-1438/ 776)  20.6 W20-W28 S18-S15
     15 1142  687 (S24W11 L=163)   195     77        0 (1110/ 654-1168/ 719)  63.4 W09-W12 S26-S21
     16  398  648 (S32E44 L=109)   836    172        1 ( 324/ 615- 480/ 695)  17.4 E51-E35 S34-S29
     17 1421  460 (S35W37 L=190)   841    141        2 (1372/ 431-1498/ 504)  -4.8 W31-W43 S37-S32
     18 1429  358 (S42W44 L=196)   815    126        2 (1388/ 325-1493/ 408)  -4.4 W38-W51 S44-S38
  

    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