Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1018 1450 (N25W04 L= 52)   576    119        1 ( 953/1421-1076/1472) -23.2 W00-W07 N22-N26
      2  332 1320 (N24E46 L=  1)  1738    322        1 ( 218/1235- 455/1423) -27.3 E63-E34 N17-N32
      3 1866  950 (S13W61 L=108)   611    110        1 (1841/ 902-1894/1004)  51.3 W58-W64 S15-S10
      4 1176  895 (S10W08 L= 55)   197     54        1 (1149/ 867-1203/ 919)  18.4 W05-W09 S11-S08
      5 1414  765 (S20W22 L= 70)   389     87        1 (1367/ 746-1454/ 790)   0.0 W18-W25 S21-S19
      6 1294  739 (S21W14 L= 62)   245     86        0 (1278/ 701-1320/ 779)  68.2 W12-W15 S23-S18
      7 1817  666 (S31W61 L=109)  1616    173        1 (1770/ 633-1888/ 717)   6.7 W54-W73 S32-S27
      8  481  501 (S27E47 L=  1)  1459    231        1 ( 428/ 395- 534/ 589)  31.0 E58-E38 S33-S21
      9  703  376 (S38E34 L= 13)   290     61        0 ( 669/ 361- 737/ 395)   8.1 E37-E30 S39-S37
  

    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