Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  775 1822 (N57E41 L=309)  1562    159        1 ( 732/1760- 845/1856)  84.3 E48-E34 N50-N61
      2  943 1599 (N39E11 L=338)   347     51        0 ( 920/1573- 972/1617)  45.0 E13-E09 N37-N40
      3  499 1495 (N28E44 L=306)   637    118        1 ( 439/1465- 549/1537) -13.0 E52-E39 N26-N30
      4 1462 1429 (N31W29 L= 19)   433    162        0 (1429/1378-1524/1484) -54.5 W27-W33 N27-N34
      5  301 1410 (N20E61 L=289)  1665    163        3 ( 250/1342- 347/1480) -28.1 E71-E53 N16-N24
      6  779 1303 (N18E18 L=332)   725    260        2 ( 674/1216- 868/1365) -25.0 E26-E11 N12-N21
      7 1808 1182 (N16W58 L= 48)   389     51        0 (1790/1147-1825/1204) -90.0 W56-W60 N13-N17
      8 1326 1080 (N07W19 L=  8)   513    134        1 (1286/1031-1372/1125)  45.0 W16-W21 N04-N10
      9 1231  695 (S17W16 L=  5)   279     82        0 (1204/ 659-1255/ 733)  68.2 W14-W16 S19-S14
     10 1488  384 (S38W46 L= 35)   375     44        0 (1462/ 368-1516/ 401)   0.0 W42-W49 S38-S36
     11 1431  359 (S40W42 L= 31)   394     39        1 (1415/ 335-1453/ 386) -45.0 W38-W44 S42-S37
     12 1079  321 (S46W12 L=  1)   363     47        0 (1052/ 302-1103/ 342) -26.6 W08-W14 S47-S44
  

    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