Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1133 1783 (N44E20 L= 17)  2643    402        5 (1013/1650-1269/1927) -22.0 E46-E04 N37-N54
      2  515 1712 (N21E59 L=338)  3269    251        8 ( 432/1623- 614/1791)  72.9 E67-E54 N14-N27
      3  876 1719 (N30E32 L=  5)   385     71        1 ( 858/1679- 894/1760) -56.3 E34-E30 N27-N33
      4  622 1713 (N23E50 L=348)   617     60        0 ( 592/1687- 651/1739)  76.0 E52-E47 N21-N25
      5  585 1514 (N10E41 L=356)  1317    384        1 ( 458/1390- 718/1659)  74.1 E45-E37 N00-N21
      6 1726 1446 (N41W39 L= 76)   542     70        1 (1698/1412-1750/1480) -53.1 W35-W41 N39-N43
      7 1677 1219 (N23W35 L= 72)   259     48        0 (1662/1184-1689/1246) -45.0 W33-W36 N21-N25
      8  353 1051 (S22E44 L=353)  1101    198        2 ( 293/ 947- 417/1102)  73.3 E50-E39 S27-S17
      9  479 1008 (S21E33 L=  4)   306     58        0 ( 444/ 991- 511/1027)  21.8 E35-E30 S23-S19
     10 1250  238 (S44W54 L= 91)   529     95        0 (1212/ 211-1304/ 271)   4.8 W47-W60 S45-S42
  

    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