Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  735 1813 (N42E42 L= 92)  1331    120        3 ( 690/1764- 786/1854) -19.7 E49-E35 N37-N45
      2  999 1684 (N35E14 L=120)   272     42        2 ( 980/1658-1020/1711) -53.1 E15-E11 N33-N37
      3  486 1642 (N26E51 L= 83)   675     83        2 ( 449/1606- 525/1674) -20.0 E57-E46 N23-N27
      4 1534 1622 (N42W29 L=163)   387     56        1 (1502/1599-1563/1643)  35.5 W25-W33 N39-N44
      5 1597 1420 (N28W33 L=167)   341     94        0 (1559/1376-1629/1455) -33.7 W29-W35 N25-N29
      6 1139 1330 (N14W02 L=136)   234     61        0 (1110/1309-1175/1346)  26.6 W00-W04 N11-N14
      7 1703 1240 (N18W42 L=176)   585    167        0 (1675/1171-1763/1292) -40.6 W39-W47 N15-N21
      8  641 1054 (S10E24 L=110)   526    195        0 ( 562/1029- 729/1079)  10.3 E28-E17 S11-S07
      9  125 1001 (S17E71 L= 63)  5394    235        7 (  76/ 887- 195/1076)  79.7 E85-E60 S23-S12
     10 1569  998 (N01W34 L=168)   292     75        0 (1533/ 976-1611/1029)  31.0 W31-W36 N00-N03
     11 1854  947 (N04W60 L=194)   528     86        1 (1826/ 904-1877/ 984) -26.6 W55-W63 N02-N06
     12  380  895 (S22E42 L= 92)   424    136        0 ( 319/ 849- 415/ 923) -26.6 E47-E38 S24-S20
     13  865  873 (S17E07 L=127)   262     50        0 ( 830/ 858- 895/ 889)  26.6 E09-E05 S18-S16
     14  168  748 (S32E71 L= 63)   769    100        1 ( 152/ 697- 183/ 800)  66.8 E73-E67 S35-S28
     15 1505  683 (S18W37 L=171)   363     76        0 (1468/ 653-1536/ 719)  53.1 W34-W38 S19-S15
     16 1544  577 (S23W44 L=178)   478     74        2 (1507/ 551-1583/ 603)  -7.1 W40-W48 S24-S22
     17  585  431 (S51E28 L=106)   786     93        1 ( 556/ 382- 609/ 480) -74.1 E30-E25 S54-S47
     18 1124  366 (S46W22 L=157)   423    103        1 (1103/ 312-1148/ 420) -48.8 W18-W26 S49-S41
  

    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