Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  676 1485 (N43E16 L=107)   390     76        0 ( 646/1459- 712/1515) -45.0 E18-E12 N40-N45
      2 1727 1413 (N13W60 L=182)   587     87        2 (1685/1385-1762/1443)  -5.7 W54-W64 N11-N14
      3 1773 1360 (N08W62 L=184)   517     60        2 (1748/1332-1790/1388)  36.9 W59-W65 N06-N09
      4  194 1127 (N26E64 L= 59)  1663    269        2 ( 128/1024- 248/1207) -32.9 E78-E56 N19-N30
      5 1860 1147 (S06W64 L=186)   895     82        2 (1824/1115-1894/1180) -35.5 W59-W68 S08-S03
      6  586  786 (N01E32 L= 91)   657    117        0 ( 558/ 728- 621/ 837) -71.6 E33-E30 S03-N03
      7 1574  650 (S28W29 L=152)   515     69        1 (1553/ 608-1596/ 691)  80.5 W27-W31 S31-S25
      8  359  614 (S07E55 L= 67)   440     43        1 ( 340/ 581- 376/ 644)  63.4 E57-E52 S09-S05
      9  397  484 (S16E59 L= 64)   468    103        0 ( 361/ 458- 416/ 533)  26.6 E64-E54 S17-S12
  

    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