Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1431 1740 (N56W32 L= 87)   625     57        1 (1391/1719-1469/1754)   4.8 W27-W39 N53-N58
      2  567 1505 (N26E39 L= 16)  1061    262        0 ( 494/1417- 630/1544) -24.8 E44-E31 N20-N28
      3 1679 1483 (N37W49 L=104)   355     62        0 (1648/1459-1705/1503)  26.6 W44-W53 N35-N39
      4  395 1448 (N20E51 L=  3)   488     90        2 ( 349/1414- 425/1474) -12.5 E56-E47 N18-N21
      5  623 1139 (N04E26 L= 28)   534    152        0 ( 552/1105- 682/1165)   0.0 E30-E21 N01-N05
      6 1698  937 (N02W46 L=101)   691     70        1 (1669/ 890-1720/ 972) -53.1 W43-W48 N00-N04
      7 1058  824 (S11W04 L= 59)   256     82        0 (1030/ 790-1090/ 859)  59.0 W02-W05 S13-S08
      8 1814  624 (S16W69 L=124)  3825    262       11 (1707/ 539-1891/ 720)  40.6 W57-W84 S21-S09
      9  803  689 (S22E11 L= 44)   320     67        0 ( 765/ 675- 848/ 706)   0.0 E13-E07 S22-S20
     10 1751  631 (S16W60 L=114)   402     66        0 (1724/ 609-1770/ 664)  63.4 W57-W62 S17-S13
     11 1487  290 (S43W57 L=112)   514     56        0 (1471/ 268-1508/ 319) -33.7 W53-W62 S44-S40
  

    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