Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1046 1876 (N57W48 L=316)   764     82        2 (1004/1847-1080/1908) -90.0 W43-W56 N53-N59
      2  759 1698 (N55W05 L=272)   281     58        0 ( 728/1684- 787/1719)   0.0 W01-W08 N53-N56
      3 1592  778 (S24W28 L=295)   685    179        1 (1509/ 728-1659/ 825)   0.0 W21-W33 S25-S22
      4 1293  681 (S21W06 L=273)   296    105        0 (1269/ 630-1332/ 716)  38.7 W03-W09 S23-S19
      5  422  677 (N00E47 L=220)   362     64        0 ( 401/ 656- 456/ 705) -45.0 E48-E44 S01-N02
      6 1359  616 (S27W09 L=276)   331     99        0 (1337/ 570-1380/ 669)  45.0 W06-W11 S29-S24
      7 1433  557 (S33W12 L=280)   445    109        0 (1379/ 541-1481/ 582) -14.0 W08-W16 S34-S31
      8 1700  503 (S49W38 L=306)  1001    134        1 (1658/ 465-1726/ 561) -48.0 W33-W45 S53-S43
      9 1080  492 (S27E12 L=255)   530    121        1 (1055/ 438-1113/ 546)  40.6 E15-E08 S29-S23
     10  498  349 (S21E60 L=207)   549     48        1 ( 473/ 325- 526/ 375)  63.4 E62-E57 S23-S19
  

    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