Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  983 1801 (N40E33 L=337)   385     62        1 ( 952/1786-1020/1814)  45.0 E35-E30 N38-N42
      2 1079 1758 (N40E23 L=346)   294     32        0 (1054/1736-1099/1777)  90.0 E25-E21 N38-N42
      3  608 1653 (N19E47 L=322)   523    116        2 ( 587/1601- 633/1706) -49.4 E51-E44 N15-N22
      4  712 1600 (N18E37 L=333)   811    211        3 ( 682/1499- 735/1694) -42.5 E44-E32 N11-N23
      5 1529 1444 (N34W19 L= 29)   900    148        4 (1502/1376-1560/1511) -71.6 W16-W22 N29-N39
      6  806 1391 (N08E23 L=347)   989    312        4 ( 701/1286- 918/1472)  68.2 E26-E18 N00-N16
      7 1225 1308 (N15W03 L= 13)   206    150        0 (1192/1255-1271/1367) -29.7 W00-W07 N13-N18
      8  278 1240 (S12E56 L=314)  1082    275        1 ( 169/1208- 386/1282)  16.9 E69-E46 S15-S07
      9 1025 1163 (N01E04 L=  6)   781    276        2 ( 933/1125-1103/1221)  29.1 E09-W00 S01-N06
     10 1628 1048 (N13W35 L= 45)   285    138        0 (1593/ 983-1652/1108) -35.5 W31-W38 N10-N15
     11  599 1063 (S15E26 L=344)   275     73        1 ( 564/1040- 632/1088)  53.1 E27-E23 S17-S13
     12 1511 1028 (N08W28 L= 37)   243     54        0 (1499/ 997-1524/1058) -56.3 W26-W28 N06-N09
     13  226  869 (S36E57 L=313)  1363    231        1 ( 142/ 791- 302/ 942)  23.2 E75-E47 S41-S29
     14  938  706 (S26W05 L= 14)  1117    332        2 ( 793/ 658-1091/ 771)  47.0 E03-W11 S32-S17
     15 1510  628 (S13W42 L= 51)   262     90        0 (1464/ 610-1553/ 643)  23.2 W38-W45 S14-S11
     16  336  566 (S54E47 L=322)   552     90        1 ( 304/ 548- 379/ 593) -11.3 E53-E40 S55-S52
     17 1365  493 (S24W39 L= 49)   308    102        0 (1312/ 485-1417/ 503)  33.7 W35-W42 S26-S22
     18 1184  303 (S42W42 L= 52)   738    114        0 (1131/ 275-1231/ 332)  74.1 W39-W45 S45-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