Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1110 1855 (N54E09 L=341)   509     69        0 (1074/1831-1138/1880)  50.2 E11-E06 N51-N58
      2  467 1708 (N33E58 L=291)   489     59        2 ( 452/1673- 480/1744) -40.6 E62-E54 N29-N35
      3 1202 1503 (N26W05 L=355)   586    193        0 (1165/1430-1228/1582) -79.7 W01-W06 N20-N31
      4 1803 1242 (N19W53 L= 42)  1100    257        1 (1763/1127-1831/1351) -77.0 W49-W55 N12-N25
      5  650  795 (S24E21 L=328)   594    205        1 ( 602/ 705- 677/ 884)  74.7 E24-E19 S29-S18
      6  469  836 (S23E35 L=315)   243     78        0 ( 434/ 812- 486/ 872)  45.0 E37-E33 S24-S20
      7 1715  652 (S17W55 L= 45)   334     47        2 (1690/ 629-1733/ 672) -11.3 W52-W57 S18-S15
      8  621  539 (S41E24 L=325)   235     71        0 ( 587/ 523- 660/ 560)  -9.5 E26-E20 S41-S39
      9 1119  185 (S62W35 L= 25)   851     48        2 (1088/ 153-1145/ 219) -31.0 W30-W42 S64-S58
  

    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