Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  883 1874 (N56E24 L=288)   951    111        2 ( 835/1833- 943/1897)  16.7 E30-E19 N51-N59
      2  905 1729 (N41E15 L=297)   798    124        1 ( 850/1694- 962/1775) -24.4 E20-E09 N38-N45
      3  944 1643 (N35E10 L=302)   264     82        0 ( 934/1599- 955/1684)  90.0 E11-E09 N31-N37
      4 1118 1572 (N31W03 L=315)   384    146        0 (1097/1501-1138/1639) -78.7 W01-W04 N25-N35
      5  376 1585 (N28E56 L=256)   416     83        0 ( 350/1550- 404/1612)  53.1 E58-E53 N25-N30
      6  506 1561 (N27E41 L=271)   301     79        0 ( 484/1522- 533/1596)  90.0 E42-E40 N24-N29
      7  633 1566 (N28E31 L=281)   451    108        0 ( 583/1525- 671/1590) -26.6 E35-E27 N24-N29
      8 1204 1374 (N18W09 L=321)   235     70        0 (1174/1358-1245/1386)  14.0 W07-W11 N16-N18
      9  677  825 (S18E21 L=291)   287     72        0 ( 659/ 784- 700/ 867)  68.2 E22-E19 S20-S15
     10  468  763 (S23E37 L=275)   933    212        1 ( 395/ 691- 563/ 816)  25.0 E43-E28 S27-S19
     11  840  733 (S23E10 L=302)   188     65        0 ( 803/ 722- 875/ 744) -11.3 E12-E07 S23-S22
     12  529  526 (S38E36 L=276)   295     69        0 ( 493/ 504- 559/ 546)  20.6 E40-E32 S39-S36
     13  807  495 (S39E13 L=299)   626     91        3 ( 768/ 464- 852/ 538)  29.7 E16-E09 S41-S36
  

    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