Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1208 1554 (N30W16 L=132)  1507    346        1 (1070/1455-1324/1649) -43.0 W07-W22 N22-N37
      2 1452 1586 (N31W36 L=152)   742    179        0 (1379/1535-1531/1635) -37.9 W30-W40 N27-N34
      3  857 1581 (N34E08 L=108)   327     51        1 ( 829/1560- 890/1603) -21.8 E10-E05 N32-N36
      4  552 1472 (N29E31 L= 86)   312    114        0 ( 504/1447- 608/1498)  20.6 E34-E26 N27-N30
      5  441 1363 (N23E38 L= 78)   992    164        1 ( 408/1306- 496/1430) -60.9 E41-E33 N18-N27
      6  921 1382 (N20E04 L=112)   224     55        0 ( 889/1366- 945/1404)   0.0 E06-E02 N19-N21
      7  183 1352 (N24E67 L= 49)  1341    106        4 ( 154/1300- 209/1400)  63.4 E72-E63 N21-N27
      8  798  820 (S13E15 L=101)   207     62        0 ( 762/ 794- 826/ 837)  26.6 E17-E13 S14-S11
      9  860  631 (S25E13 L=103)   556    158        0 ( 840/ 567- 905/ 710) -84.3 E14-E11 S29-S19
     10 1640  545 (S35W46 L=162)  2112    344        4 (1515/ 443-1747/ 667)  47.0 W36-W54 S42-S27
     11 1403  473 (S39W26 L=142)   716    118        2 (1366/ 414-1433/ 525) -48.8 W22-W30 S43-S35
  

    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