Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  872 1616 (N46W10 L= 72)   347     47        0 ( 853/1585- 896/1640)  31.0 W07-W12 N44-N48
      2  951 1611 (N43W15 L= 77)   306     50        2 ( 920/1587- 976/1630) -45.0 W13-W17 N41-N45
      3 1569 1474 (N16W49 L=111)   341    102        1 (1556/1429-1591/1525)  39.8 W45-W51 N13-N18
      4 1241 1464 (N25W27 L= 89)   222     86        0 (1190/1448-1274/1481) -18.4 W23-W29 N24-N26
      5 1448 1368 (N14W35 L= 97)   509    125        3 (1384/1344-1507/1388) -20.6 W31-W39 N12-N15
      6  779 1199 (N23E10 L= 52)   198     55        0 ( 757/1173- 800/1225)  33.7 E11-E08 N21-N24
      7 1747 1099 (S10W48 L=110)   554    112        1 (1707/1049-1785/1141)  20.6 W43-W52 S11-S08
      8 1659 1034 (S11W39 L=101)   361     93        1 (1632/ 995-1686/1080) -80.5 W37-W40 S13-S07
      9  468  984 (N18E36 L= 27)   585    168        0 ( 415/ 935- 535/1043)  15.3 E40-E29 N15-N19
     10  257  935 (N19E56 L=  6)   436     84        1 ( 228/ 906- 291/ 976) -39.8 E58-E52 N16-N21
     11  776  930 (N08E16 L= 46)   269     92        1 ( 760/ 885- 795/ 967)  68.2 E18-E15 N05-N10
     12  155  877 (N17E71 L=351)   769     89        2 ( 137/ 834- 174/ 921) -71.6 E73-E68 N13-N19
     13 1316  836 (S11W12 L= 74)   402    125        0 (1292/ 792-1341/ 887)  71.6 W09-W13 S13-S07
     14 1598  779 (S23W29 L= 91)   572    144        1 (1553/ 718-1640/ 833)  26.6 W24-W33 S25-S21
     15 1686  738 (S29W36 L= 98)   351     61        0 (1664/ 710-1706/ 774)  45.0 W33-W37 S30-S27
     16 1169  656 (S18E01 L= 61)  1090    381        5 (1007/ 622-1358/ 688) -21.8 E10-W10 S21-S13
     17 1710  590 (S41W38 L=100)   518     48        2 (1681/ 567-1737/ 619) -20.6 W34-W42 S43-S38
     18  403  435 (S14E63 L=359)   968    138        1 ( 382/ 374- 419/ 505)  42.0 E68-E57 S18-S09
     19  687  412 (S21E40 L= 22)   361     83        1 ( 669/ 377- 709/ 453)  71.6 E41-E38 S23-S17
  

    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