Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1109 1763 (N47W38 L=117)   329     72        0 (1083/1748-1145/1784) -63.4 W35-W40 N45-N49
      2  902 1690 (N50W17 L= 96)   428    126        1 ( 875/1633- 949/1722)  20.0 W11-W22 N46-N51
      3 1667 1628 (N19W70 L=150)  1722    174        3 (1625/1554-1705/1687)  48.8 W65-W76 N14-N22
      4 1043 1431 (N29W13 L= 93)   357    148        0 ( 981/1390-1098/1469) -58.0 W10-W15 N24-N32
      5  487 1358 (N40E28 L= 52)   584    152        1 ( 455/1291- 521/1429) -70.0 E31-E25 N34-N45
      6  798 1270 (N26E07 L= 73)   540    173        3 ( 707/1244- 874/1294) -31.0 E12-E02 N23-N29
      7  419 1199 (N31E36 L= 44)   263     78        1 ( 390/1167- 436/1228) -76.0 E38-E34 N29-N33
      8 1893 1168 (S12W66 L=146)   946     82        0 (1847/1148-1935/1186)  -7.1 W59-W75 S14-S10
      9  459 1140 (N27E33 L= 47)   451     96        5 ( 408/1109- 500/1173) -14.0 E37-E29 N24-N28
     10  553 1078 (N21E27 L= 53)   394     93        1 ( 532/1044- 578/1114) -76.0 E28-E25 N19-N23
     11  355  953 (N18E45 L= 35)   828    289        0 ( 274/ 853- 471/1032) -45.0 E51-E37 N10-N24
     12  626  943 (N12E25 L= 55)   258     96        0 ( 594/ 902- 664/ 983) -71.6 E26-E23 N08-N14
     13  501  874 (N11E35 L= 45)   770    136        4 ( 473/ 798- 541/ 932)  63.4 E38-E32 N06-N14
     14 1545  770 (S23W24 L=104)  1011    388        0 (1408/ 669-1654/ 912)  17.7 W13-W35 S25-S17
     15 1380  807 (S15W14 L= 94)   236     46        2 (1357/ 785-1410/ 824)  18.4 W13-W16 S16-S13
     16 1799  758 (S33W48 L=127)   829    107        1 (1751/ 719-1840/ 792)   5.2 W41-W53 S34-S31
     17 1230  620 (S22W01 L= 81)   453    196        0 (1172/ 562-1286/ 706)  31.0 E04-W06 S24-S18
     18  725  370 (S24E40 L= 40)  1411    348        1 ( 619/ 251- 809/ 465)  17.8 E58-E30 S31-S18
  

    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