Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1155 1739 (N44W05 L= 63)   283     41        0 (1132/1719-1181/1755)  26.6 W02-W06 N42-N45
      2 1352 1669 (N40W20 L= 79)   499     60        0 (1329/1640-1385/1692)  31.0 W18-W23 N37-N42
      3  909 1621 (N33E13 L= 46)   788    153        0 ( 855/1551- 960/1671) -45.0 E16-E08 N27-N36
      4 1453 1626 (N38W28 L= 86)   549     71        0 (1431/1580-1478/1656)  63.4 W26-W30 N34-N40
      5 1672 1507 (N31W45 L=103)  1808    289        2 (1569/1438-1791/1560) -11.9 W35-W55 N27-N34
      6  790 1307 (N11E16 L= 42)   239     35        0 ( 759/1291- 814/1325)   0.0 E18-E14 N09-N12
      7  219 1260 (N07E59 L=359)   387     77        0 ( 190/1234- 241/1282) -31.0 E63-E56 N05-N08
      8  445 1196 (N03E38 L= 20)   231     54        0 ( 413/1175- 475/1220)  36.9 E40-E35 N01-N04
      9  840  689 (S26E10 L= 48)   275     78        0 ( 815/ 657- 867/ 735)  59.0 E11-E08 S28-S23
     10  401  572 (S36E45 L= 13)   884    170        1 ( 338/ 543- 478/ 610) -20.0 E50-E38 S37-S32
  

    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