Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  734 1303 (N24E18 L=195)   252     80        0 ( 711/1269- 755/1335) -63.4 E19-E16 N21-N25
      2  725  935 (N01E19 L=194)   239     63        0 ( 702/ 919- 754/ 958)  18.4 E20-E17 N00-N02
      3 1245  893 (S05W13 L=226)   601    201        2 (1154/ 860-1324/ 933) -26.6 W07-W17 S07-S01
      4 1561  830 (S11W34 L=248)   361     78        1 (1519/ 816-1607/ 850)  -9.5 W31-W37 S12-S09
      5  207  778 (S09E65 L=148)  1078    101        2 ( 182/ 736- 237/ 827)  71.6 E68-E60 S11-S05
      6 1669  703 (S21W45 L=258)   528    110        1 (1637/ 652-1696/ 749)  56.3 W42-W47 S23-S17
      7  652  608 (S20E28 L=185)   671    178        1 ( 608/ 547- 718/ 662)  26.6 E32-E22 S23-S16
      8 1326  542 (S28W18 L=231)   283     82        1 (1307/ 500-1341/ 579) -80.5 W16-W19 S31-S25
      9  574  508 (S27E36 L=177)   627     90        2 ( 549/ 471- 616/ 545) -45.0 E39-E33 S29-S23
     10  561  376 (S37E44 L=169)   731    124        1 ( 546/ 310- 591/ 424)  51.3 E50-E38 S42-S32
  

    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