Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1536 1728 (N30W62 L= 72)   983    115        1 (1496/1682-1582/1762)   0.0 W55-W71 N27-N31
      2 1792 1247 (S03W56 L= 65)   491     57        1 (1764/1222-1821/1275)  21.8 W52-W58 S04-S01
      3 1905 1186 (S09W67 L= 77)   537     47        4 (1890/1161-1924/1210)  18.4 W64-W70 S10-S07
      4  206 1146 (N27E57 L=313)   399     53        0 ( 178/1126- 238/1166)   8.1 E60-E52 N26-N28
      5 1313 1018 (S04W16 L= 26)  1972    572        5 (1124/ 837-1448/1151) -53.6 W06-W20 S16-N05
      6 1810 1124 (S10W54 L= 64)   379     50        0 (1788/1099-1841/1142)  -9.5 W51-W57 S11-S08
      7 1587 1083 (S07W35 L= 44)   609    143        1 (1545/1015-1620/1135) -82.9 W32-W36 S10-S02
      8  413 1114 (N21E36 L=333)   257     62        0 ( 380/1093- 444/1132) -36.9 E38-E33 N19-N22
      9 1669 1011 (S13W39 L= 49)   930    290        3 (1593/ 935-1726/1094)  21.0 W32-W45 S15-S09
     10  752 1036 (N10E15 L=355)   629    166        1 ( 662/1008- 824/1062) -29.1 E20-E11 N07-N12
     11 1680  821 (S24W37 L= 47)   524    105        0 (1660/ 778-1706/ 881) -90.0 W35-W39 S26-S20
     12 1772  791 (S29W46 L= 56)   905    148        0 (1706/ 769-1841/ 824)  -4.1 W39-W53 S30-S27
     13 1549  735 (S26W26 L= 36)   451    111        1 (1534/ 681-1565/ 800) -90.0 W24-W26 S29-S21
     14  905  716 (S11E14 L=356)   646    150        3 ( 857/ 685- 968/ 754) -45.0 E16-E10 S13-S08
     15  764  475 (S23E30 L=340)  3410    815        6 ( 513/ 247-1006/ 689)  16.9 E65-E09 S33-S14
     16 1556  555 (S38W25 L= 35)  2735    272        7 (1475/ 447-1638/ 662)  52.4 W18-W32 S46-S32
     17  477  582 (S10E45 L=325)   421    150        0 ( 409/ 563- 535/ 619)  -5.2 E51-E39 S10-S08
     18 1415  478 (S39W12 L= 22)   294     54        1 (1387/ 456-1439/ 501)  26.6 W10-W14 S40-S37
  

    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